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Abstract 

This paper presents a detailed assessment of the ability of the 240 Skyrme interaction parameter 
sets in the literature to satisfy a series of criteria derived from macroscopic properties of nuclear 
matter in the vicinity of nuclear saturation density at zero temperature and their density depen- 
dence, derived by the liquid drop model, experiments with giant resonances and heavy-ion collisions. 
The objective is to identify those parameterizations which best satisfy the current understanding 
of the physics of nuclear matter over a wide range of applications. Out of the 240 models, only 
16 are shown to satisfy all these constraints. Additional, more microscopic, constraints on density 
dependence of the neutron and proton effective mass beta-equilibrium matter, Landau parameters 
of symmetric and pure neutron nuclear matter, and observational data on high- and low-mass 
cold neutron stars further reduce this number to 5, a very small group of recommended Skyrme 
parameterizations to be used in future applications of the Skyrme interaction of nuclear matter 
related observables. Pull information on partial fulfillment of individual constraints by all Skyrme 
models considered is given. The results are discussed in terms of the physical interpretation of the 
Skyrme interaction and the validity of its use in mean-field models. Future work on application of 
the Skyrme forces, selected on the basis of variables of nuclear matter, in Hartree-Fock calculation 
of properties of finite nuclei, is outlined. 

PACS numbers: 21.30.Fe, 21.65.Cd, 21.65.Ef, 26.60.Kp 
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I. INTRODUCTION 



Empirical properties of infinite nuclear matter can be calculated using many different 
theoretical approaches. The most microscopic ones start from a realistic two-body free 
nucleon-nucleon (NN) interaction with parameters fitted to NN scattering phaseshifts in 
different partial wave channels and to properties of the deuteron [1'] . Taking these bare NN 
interactions as input into many-body formalism, such as the relativistic Dirac-Bruckner- 
Hartree-Fock (DBHF) approximation and its non-relativistic counterpart BHF , varia- 



tional methods 
(SCGF) 
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,|5|, Correlated Basis Function jo] ^elf-Consistent Greens Function models 
Quantum Monte Carlo techniques 9l-|l4| and Chiral Effective Field Theory 
an effective NN interaction, which includes the effect of the medium, is derived and 
the many-body problem approximately solved. 

Various many-body approaches typically lead to an over prediction of the saturation den- 
sity po of symmetric nuclear matter (SNM), at which the binding energy per nucleon reaches 
its maximum, and of the corresponding maximum binding energy Eo{p=Po) [l7]- There are 
many ways of estimating the experimental value of po, including different variants of the 
liquid drop models, optical model of NN scattering, muonic atoms, and Hartree-Fock (HF) 
calculation of nuclear density distributions (see e.g. 18| and references therein). The range 
of results is rather broad but a consensus value is po=0.17±0.03 fm^'^. The empirical value 
Eo per nucleon ~ 16 MeV can be extracted from the semi-empirical mass formula or from the 
extrapolation of binding energies of heavy nuclei. Theoretical calculation of saturation prop- 
erties of SNM is not only dependent on the choice of the bare NN interaction but also on the 
method of treatment of many-body effects. For example, if the BHF approximation is used. 



Eo and po are correlated within a narrow band [17|, |19| . Two main approaches have been 
suggested to improve the theoretical calculation of saturation properties of SNM, the most 
frequently used being the inclusion of 3-body (NNN) forces. As the form of these forces is 
unknown, different ad hoc parameterizations have been used, dependent on additional vari- 
able parameters that need to be fitted to account for the delicate balance between the strong 
(NNN) attraction and (NN) repulsion at short distances. Alternatively, DBHF calculations 

20|. Another possi- 



have been shown to be effective without the need for NNN forces 
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)ility is to treat the scalar and vector densities in the Walecka relativistic mean-field model 



2l| as equal 22|]. However, the systematic deviation of all theoretical predictions from the 
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expected empirical values of Eo and po remains a problem. An interesting suggestion by 
Dewulf et al. ^\ implies that treatment of short-range correlations in nuclear matter in the 
SCGF model brings the saturation density closer to the empirical value than do current 
BHF calculations. Careful examination of the effect of long-range correlations on nuclear 
saturation properties through coherent pion-exchange contributions to the binding energy 
of nuclear matter is equally important. Such correlations are not present in finite nuclei and 
a question arises as to how applicable are liquid drop model predictions, based on properties 
of finite nuclei, to saturation properties of nuclear matter. 

Quite a different perspective on treating many-nucleon systems is to use effective den- 
sity dependent NN and NNN interactions instead of realistic ones. The pioneering models of 



Kohler 



23| , Brink and Boeker 2J] , Moszkowski 25| , Skyrme 



therin and Brink (see 



26 



27( 1 and references therein) and Gogny 



aarther developed by Vau- 
initiated this approach. 



widely used today. The basic idea is to parameterize the NN and NNN interactions by 
zero range (Skyrme model), short finite range (Gogny model) and indefinite range (Separa- 



ble Monopole model (SMO) [29 



30|) density dependent functionals to describe the ground 



state properties of finite nuclei and nuclear matter. In this scenario, the microscopic de- 
tails of NN and NNN forces, such as meson exchange, are not explicitly considered and 
all the physically relevant information is carried by the parameters of the density depen- 
dent phenomenological forces which include the spin, orbital angular momentum and isospin 
couplings. The drawback of this approach is that the parameterization of such forces is not 
unique and there exist, in principle, an infinite number of parameter sets, fitted to ground 
state properties of (doubly- or semi-magic) stable nuclei, fission barriers, energies of giant 
resonances and symmetric and asymmetric nuclear matter (ANM). This situation arises in 
part because there is no unambiguous connection between individual parameters, or groups 
of parameters, of these forces with particular physical properties of the many-body nuclear 
system. Many parameters are strongly correlated. 

It is obviously desirable to constrain the parameterizations of effective density dependent 
forces as much as possible. The strategy chosen in this work is to concentrate first on 
application of the forces to modelling different variants of infinite nuclear matter. Although 
nuclear matter is an idealized medium, and all its properties, derived from experiments 
indirectly in a model dependent way, are empirical quantities, it offers an important insight 
into specific parts of the phenomenological interaction and has important applications in the 
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theory of heavy-ion colhsions (HIC) and the physics of neutron stars. This paper focuses on 
the Skyrme interaction and is the first part of a series in which nuclear matter constraints 
will be applied to each class of effective density dependent interaction. It will be followed 
by investigation of Gogny and SMO forces and interactions used in relativistic mean-field 
(RMF) models. After this work is complete, the implications of the consequences for the 
theory of finite nuclei will be investigated. 

We use the most up-to-date constraints on nuclear matter properties which go much 
beyond the minimal conditions on SNM (saturation density, binding energy, incompressibil- 
ity and symmetry energy at saturation). New data from HIC, giant monopole and dipole 
resonance experiments as well as new observational data on neutron stars provide new con- 
straints on the performance of individual Skyrme parameterizations in nuclear matter. 

The set of eleven macroscopic constraints used in this work have been mainly derived 



from experimental data, on the assumption of the validity of the liquid-drop model [3l|, |32| . 
and concern properties of SNM at and close to the saturation point. Studies of dilute Fermi 
gas provide constraints on low density pure neutron matter (PNM) equation of state (EoS). 
Dynamical models of HIC further constrain the density dependence of pressure in SNM 
and PNM at subsaturation density and extrapolate these constraints to higher densities. 
Mean- field Hartree-Fock -|- Random Phase Approximation (RPA), used in calculating giant 
resonance excitation energies both in relativistic and non-relativistic models, provide a final 
group of constraints on the incompressibility of nuclear matter and its density and symmetry 
dependence. 

In addition to the above macroscopic constraints, several more microscopic constraints 
are employed. These include the density dependence of the nucleon effective masses in beta- 
equilibrium matter (BEM), the Landau parameters for SNM and PNM and of the symmetry 
energy. Observational data on cold non-rotational high-mass and low-mass neutron stars 
provide a final group of constraints. 

In this work we consider 240 Skyrme parameter sets, currently available in the literature, 
and critically compare their predictions for a wide variety of properties of SNM, PNM, asym- 
metric matter with fixed ratio proton fraction (ANM) and BEM with all available constraints 
in the density range from ~0.1po to 3po5 estimated on the best available experimental and 
theoretical grounds. The range of applicability of the Skyrme force is a very important issue 
which is often mishandled. 
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The Skyrme interaction, originally constructed for finite nuclei and nuclear matter at 
saturation density, is a low momentum expansion of the effective two-body NN interaction 
in momentum space, and both the lower and upper limits of its validity are not firmly 
established. The important point about the Skyrme interaction is that some correlation 
effects are included through its parameters. Thus, although formulated as zero-range in 
coordinate space [27|, it exhibits some finite-range features [33 1. 

For finite nuclei, the best evidence for a lower limit derives from the fact that Skyrme 
models reasonably predict the observed abrupt decrease of density at the nuclear surface, 
and neutron and proton mean-square radii. The sensitivity extends down to about O.lpo- 

In uniform SNM the Skyrme interaction has been used to make prediction of the appear- 
ance of light clusters (deuterons, tritons, ^He and alpha-particles in hot matter in the region 
of density of 0.6 - 1.25po 3J]). More recently, abundance of light clusters with A <13 in 
supernova envelopes at finite temperature was calculated at density range 0.01 - O.Spo using 



the Skyrme functional 



35 . T 



le "pasta" phase, predicted in neutron star and supernova 



matter in a variety of models 



36 



37| . was successfully modeled in supernova matter in a 



density range 0.25 - 0.75po in self-consistent HF+BCS calculation with SkM* and SLy4 
Skyrme interaction [38]. 

PNM has been mainly studied as an approximation to a low density Fermi gas. Schwenk 
and Pethick 



39| explored, in a model independent way, the neutron matter EoS at densities 
0.0125 - 0.125po, and Epelbaum et al. |40| calculated the ground state energy of dilute 
neutron matter at next-to-leading-order in lattice chiral effective field theory in the density 
range 0.02 - O.lpo- Quantum Monte Carlo techniques have been applied to low density PNM, 



41 



42| . To our knowledge. 



providing a constraint for on the EoS up to saturation density 
there has not yet been a detailed study of the applicability of the Skyrme interaction at 
these low densities in PNM. Such a study is of a particular interest as it may be one of the 
best ways to model the crust of neutron stars. 

The upper density limit of validity of the Skyrme interaction reflects the fact that at 
higher densities relativistic effects should be increasingly important. The appearance of 
heavy strange baryons and mesons in the matter is ultimately inevitable. Due to Pauli 
blocking, the chemical potential of the neutrons increases rapidly with density. At some 
point, it becomes energetically favourable for the system to let the neutrons undergo a 
strangeness changing weak decay, which replaces them by hyperons, for which the Fermi sea 



6 



is not yet filled. From the difference of mass between the neutron and its strange partners it 
follows that the critical density, at which hyperons should appear, is 2-3po- Using a nucleon- 
only Skyrme interaction beyond this density can be expected to yield misleading results. 
This is discussed later in connection with high-mass neutron star models. 

Taking all the above pieces of evidence into account, we adopt O.Olpo < P < 3po as the 
range of validity of the Skyrme interactions considered in this work. 

The paper is organised as follows. A brief description of the Skyrme interaction, to- 
gether with definition of the variables used in this work, is given in Sec. [Tll Classification of 
the macroscopic constraints and discussion of their origin and applicability range forms the 
content of Sec. Illli Sec. IIVI presents a comparison of predictions of those Skyrme parameteri- 
zations which satisfy the macroscopic constraints with further microscopic and observational 
constraints. The results are discussed and summarized in Sec. |V] and conclusions are pre- 
sented in Sec. |Vll 



II. SKYRME MODELS 



Since the original work by Skyrme in the fifties [26| and the Vautherin and Brink 27 1 
parameterization of the original interaction in early seventies, considerable effort has been 
invested in the application of this density dependent effective interaction both to ground 



state properties of finite nuclei and to nuc 



ear matter in the framework of the mean-field 



Hartree-Fock approximation (see e.g. 43, 



44l | for recent reviews). The advantage of the 



structure of the Skyrme density functional is that it allows analytical expression of all vari- 
ables characterising infinite nuclear matter 



45 



47| . Such structure can also be constructed 



from non-relativistic versions of the relativistic point-coupling models |48l-l50| . In the fol- 
lowing, we introduce the various physical quantities and give expression for each in terms of 
the Skyrme parameters. The general expression for the energy per particle of infinite ANM, 
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defined in terms of tlie energy density 8 and particle number density p, is given 
3R= /3ir'^"' 



p lOM 



witfi 



p^/^i^5/3 + |p[2(xo + 2)-(2xo + l)i/: 



2J 

2\ 2/3 



+ ^ E^3.P'^''''[2(^3. + 2) - (2x3. + l)H,] + ^ ( — 1 P'" («^5/3 + hH,„) 

i=l ^ ^ 



+ 



1 




40 


(1 






40 





2\ 2/3 



2\ 2/3 



-,5/3+(5 



,5/3+7 



t4(a;4 + 2)if5/3 - ^4(2:4 + 2)-^8/3 



t5(x5 + 2)i75/3 + t5(a;5 + -)i^8/3 



(1) 



a = ti(xi + 2) +t2(a;2 + 2), 

h = ^ [t2(2x2 + 1) - ti(2xi + 1)] , and 



(2) 
(3) 
(4) 



q. (II]) includes tfie summation over index i in tlie 



wfiere y = Z/A is the proton fraction, 
third term introduced by Agrawal et al. 5l| and additional terms involving ^4, X4, and ts, X5, 
used by Chamel et al. [5^. The great majority of the parameterizations referred to in this 
work do not include these terms. Parameterization without (with) these additional terms 
are regarded as "standard" ("non-standard") in this paper. We note that there are several 
other parameter sets which parameterize the density dependence of the Skyrme functional 
in non-standard ways different from those considered here. These forces have been 

reported to have problems at higher density nuclear matter 5l|] and have not been included 
in the present study. 

All quantities referred to in this work have been obtained based on Eq. ([1]) and are given 
below. 

Eq. (Il]) leads to an in-medium effective nucleon mass M* in ANM 

M* = M |i75/3 + \^P [{a + U{xi + 2)p^ + U{x^ + 2)p^) H,,^ 



+ (6-t4(x4 + ^)p^ + t5(a;5 + ^)p^)i^8/3 



-1 



(5) 



with M being the free nucleon mass. 
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The pressure, defined as P—p^^^^^, is given as 



5M 



p'"H,/, + |p^[2(xo + 2) - (2x0 + l)H, 



1 3 1 /Stt^X^^^ 

+ 4^ E + 1)P''^^'[2(^3. + 2) - (2x3. + + 81^) P'^' («^5/3 + hHs/^) 



1 /37r^ 
+ 40 



2\ 2/3 



(5 + 2,5)p 



+ 



1 /37r^ 



40 



2\ 2/3 



(5 + 37)p= 



^4(2^4 + 2)-ff5/3 - ^4(2^4 + 2)^8/3 



t5{x5 + 2)H5/3 + t5{x5 + -)Hs/3 



(6) 



The volume incompressibihty of ANM at saturation density is calculated as derivative of 

pressure with respect to number density p 

'd^S/p' 



K = 



dp"^ 



-(f) 

= ^ y—j p'/'^5/3 + ^p[2(xo + 2) - (2x0 + IW 

3 3 

+ T7:Y.^M{(Ji + l)(<Ji + 2)p'^^+^[2(x3, + 2) - (2x3i + l)//2] 



1=1 



+ 3 



37r^ 



2\ 2/3 



p'''\aH^,^ + hH^,^) 



+ 



3 /37r^ 



40 



2\ 2/3 



3 /37r^ 
+ 40 



2\ 2/3 



(5 + 35)(8 + 2,5)p-^+'[U{x^ + 2)i/5/3 - U{x^ + -)i/8/3] 
(5 + 37)(8 + 37)pi+^[i5(x5 + 2)H^/^ + ^5(^5 + ^)i^8/3]. 



(7) 



Finally, the third derivative of the energy per particle in ANM at saturation density, also 
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called the skewness coefficient, is expressed as 

9 ^ 

+ 5^^3it^i(fXi + l)icTi - l)p'^»+i[2(x3i + 2) - (2x3i + l)H2] 
lb ^ — ^ 

1=1 

3 /Svr^^'/^ 



4 V 2 



p5/3(aff5/3 + &^^8/3) 
2\ 2/3 



+ ^ + ^^^^^ + ^^^^^^ " l)P^'"'[i4(x4 + 2)i/5/3 - t4(x4 + ^)i^8/3 

+ ^(?) ^ (2 + 37)(5 + 37)(37-l)pi+^[t5(x5 + 2)i75/3 + t5(x5 + ^)i^8/: 



(8) 



In a special case for the SNM {y = 1/2), the expressions (P - E]) simplify in that 
= 1 and we obtain energy per particle EsNuip), the incompressibility -K'snm(p) and 
skewness (5snm(p)- When calculating properties of symmetric matter at the saturation den- 
sity -Esnm(Po) = Eo, the incompressibility -R'snm(Po) = Ko and the skewness (5snm(Po) = Qo, 
the second and third derivatives of the energy density with respect to number density are 
taken at p = po. Obviously, the first derivative, the pressure, is equal to zero at po- 

One of the key properties of nuclear matter is the symmetry energy, particularly important 
in modelling nuclear matter and finite nuclei because it probes the isospin part of the Skyrme 
interaction. It is defined as 

8 Qy2 \p,y=l/2 

/3^2x2/3 , 3 



5(P) 



6M 



(1^) P'^' - |(2xo + 1)P - ^ E ^3^2X3. + 1)P-^^ 
^ ^ 1=1 

24 (—J + -8(~J ''''' 

t5(5x5 + 4)pt+^. (9) 



1 /37r2N^/^ 



24 

In SNM, it is customary to define four quantities, J=5(po), L, -K'sym and Qsym related to the 
symmetry energy and its derivatives evaluated at the saturation density po- L, the slope of 
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iS, is given by 



dp J 



P=Po 

2 /q„2\2/3 oy- i 3 

^,/3 _ ^ ^^^^ _ 1 ^,3^(2X3. + + l)p-+^ 



3M 



5 A3.r^^^^ . <>vV3_ir^V'%....U...J-^ 



+ ^ J (a + 46)p^/^ - o ( — ) (5 + 35)t4X4Po^ 

2\ 2/3 



1 /3 

+ ^ j (5 + 37)t5(5x5 + ^ft\ (10) 

The curvature of the symmetry energy S at saturation density in SNM is sometimes called 
i^sym, the symmetry incompressibility. It should not be confused with K^, which is the 
isospin incompressibility, defined in Eqs. (1^ and iiTsym is given by 



K = 9o 



sym -ro ^ 



^ ^ P=Po 



3M 



^ ^ i=l 



+ ^ (^) (a + 46)P^/^ - i (^) (5 + 35) (2 + 



1 /37r2 



2\ 2/3 



(5 + 37)(2 + 37)^5(5x5 + A)pr . (11) 



24 

Finally, Qsym, the third derivative of the symmetry energy, is 

Qsyvn 27Pq 



3M 



9p- 

/q 2\^/'^ q 

( V ) '^o^' - E^3,(2X3, + l)(a, + l)a,(a, - l)p- 
^ ' 1=1 



P=Po 



+1 



5 /37r2^2/^ ^ /37r2\^/^ 



12 



a + 4&)p^/3 _ 1 ( r_ ) (5 + (2 + 35) (35 - l)Uxipi 



1+5 



1 /o 2\2/3 g 

+ ^ (^— J (5 + 37)(2 + 37)(37 - lM5x, + 4)pt\ (12) 

Using the above expressions, the density dependence of the symmetry energy can be 
expanded as a function of a; = (p — po)/3po, 

S = J + LX + ^KsymX^ + ^QsymX^ + 0{x^) . (13) 
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Similarly, in SNM, the density dependence of the energy per particle E, Eq. ([T]), is 
sometimes expressed as an expansion in a Taylor series around po 

EsMp) = Eo + Ikx^ + IQox' + Oix^). (14) 
2 o 



where Eq = Esnm(Po) is the energy per particle at the saturation density po and 

5sn: 
dp 



K = 9pl{^^^^] (15) 



P=Po 



and 



In ANM with asymmetry /3 = (A^ — Z)/A = (1 — 2y), the energy per particle E can be 
expanded around a new, isospin dependent, saturation density Po(/3) ~ Po(l — '^i.L/ Ko)[5'^) 



57|: 



(17) 

where the expansion coefficients are given as 

EM(3)) = E, + J(3^ + 0{(3^), (18) 

K^iPom = K,+ (^i^sym - 6L - I^L^ f3^ + 0{/3^), and (19) 

QoiPom = Qo+ (Qsym " 9l|^^ (3^ + 0{(3'). (20) 

The coefficient of the second term in Eq. (fT9l) 

i^,,, = (K^ym - 6L - (21) 



determines the isospin dependence of incompressibility at saturation density Po{f3). Strictly 
it is the volume part Kr^v of the isospin incompressibility Kr Eq. fl24l) . which plays an 
important role in analysis of data from giant monopole resonance. It does not include 
surface effects, as discussed in Sees. IIII Al and IIII CI 

III. MACROSCOPIC CONSTRAINTS 

It is important to keep in mind that different Skyrme parameterizations were often con- 
structed with emphasis on a certain selection of data on finite nuclei. For example, the 
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BSk family were fitted to experimental nuclear masses, SkM* to binding energies of finite 
nuclei and actinide fission barriers, the Ski family to isotope shifts in the Pb region, and 
the SLy family to properties of neutron matter, neutron stars and the ground state vari- 
ables of neutron-heavy nuclei. Although all Skyrme forces are usually fitted to reproduce 
well the saturation energy and density of symmetric nuclear matter, they differ significantly 
in other characteristics of symmetric and pure neutron matter, in particular their density 
dependence. 

We examine in this section eleven constraints on properties of nuclear matter, out of 
which four are related to SNM, two to PNM, and five involving both SNM and PNM. The 
constraints are listed in Table HI 



A. Symmetric nuclear matter 

Infinite nuclear matter, composed of the same number of protons and neutrons without 
Coulomb interaction, does not exist in nature. Nevertheless, it has become an important 
theoretical laboratory for the investigation of physical quantities relevant for the modelling 
of heavy nuclei and nuclear matter in astrophysical compact objects. As stated above, the 
saturation density Po and the binding energy per nucleon Eq = £/po are reasonably well 
established. In this work we focus on two other physical quantities of SNM at saturation 
density: The incompressibility and the skewness coefficient Qo- We note that some 
authors use K' = —Qo as the skewness coefficient (see Eq. ([8])). 

The determination of Ko, and other related parameters of nuclear matter from experi- 
mental data on Giant Monopole Resonance (GMR) on finite nuclei, has been a long-standing 



problem (see e.g 58l-l60|) which has not been fully resolved to this day 61|. 

There are basically two procedures that can be used to extract the information. One is 
to analyse the A dependence of the compression modulus of a finite nucleus Kj.^ 

Kj, = (M/h') <R'> El^^ (22) 

obtained, using sum-rule arguments, from the measured energy of Giant Monopole Reso- 



nance -Egmr in spherical nuclei 



M is the nucleon mass and R is the rms matter radius 



of the nucleus with mass number A. Using the leptodermous expansion of the energy per 
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particle in the droplet model (the mass formula), an expression for can be found 



+ K,/32 + ir,oui^ + --- , (23) 



where the isospin incompressibility K^. consists of two components 58l. |62|. |63| 



Kr = Kr^, + (24) 

where Kr^^, (-^r,s ) determine the volume (surface) isospin incompressibility. K^o\, -ft'surf, 
J^curv, Kj. and -ft'coui are second derivatives of the coefficients of the volume, surface, curva- 
ture, isospin, and Coulomb terms with respect to the radial coordinate of the nucleus in the 
mass formula. If the K^o\ term is identified as the incompressibility of infinite nuclear matter 
then a link between Kp^ and can be used to determine K^- However, it has been ar- 
gued that the macroscopic analysis suffers from potentially serious drawbacks, including the 
uncertainty of the validity of the i^voi ~ Ko assumption, the weak A dependence, limitations 
of Eq. fl23l) to small oscillations, questions about convergence of the expansion Eq. fl23l) . and 
consideration of possible anharmonicities of the breathing mode, especially for light nuclei 



60|. 



Another route is to rely on microscopic calculations within a Hartree-Fock mean-field 
approximation for static properties and RPA for excitations, and to use the same model 
framework for calculation of the SNM. Such an approach allows a consistent determination 
of the breathing mode energy -Egmr and parameters of nuclear matter within the same 
framework. However, the microscopic approach has the disadvantage that the results are 
dependent on the choice of the effective interaction used in the models. The macroscopic 
approach is based only on the assumption that the liquid drop model of the nucleus is valid 
and that the leptodermous expansion fl23l) converges reasonably fast. 

The value of most frequently used today is based on microscopic analysis of GMR 



64l | used measured EO strength distribution in ^"^Ca, ^°Zr, ^^®Sn, 



data. Youngblood et al. 
^^^Sm and ^osp]^ g^j^^ Gogny interaction based on calculations by Blaizot et al. i60], which 
took into account pairing and anharmonicity in lighter nuclei. The deduced value of 



231 ±5 MeV. This result was in good agreement with the value obtained by Farine et al. 65 1 



Ko = 240 MeV with generalized Skyrme forces. Myers and Swi<itecki 66| used a model based 



on a semiclassical Thomas-Fermi approximation with a short range Yukawa effective force 
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between fermions. They obtained the value Ko= 234 MeV by fitting nuclear binding energies 



and diffuseness to experimental data. Later, Colo and co-workers [67|] argued that there is 
not a unique relationship between the value of Ko associated with an effective force and the 
GMR energy predicted by that force in non-relativistic models. They built a new class of 
Skyrme forces and found that 230< Ko <250 MeV. Agrawal et al. addressed the well-known 
discrepancy between predictions of the value of Ko in non-relativistic and relativistic models 
(see [68j] and references therein) and concluded that the discrepancy, thought to be ~ 20%, 
can be much smaller for an appropriate choice of effective interactions. Todd-Rutel and 

n 

Piekarewicz |69J| , using a new relativistic model FSUGold with two new parameters, causing 
softening both of the EoS of SNM and the density dependence of the symmetry energy, 
predicted Ko= 230 MeV in RMF calculations. However, new data on GMR in Sn isotopes 
7ol . 7l| re-opened the question of Ko, as models which successfully predicted GMR energies 



in ^°Zr, ^^^Sm and ^°^Pb could not reproduce GMR energies reported for ^^'^Sn. We note 
that the new GMR energies for ii2,ii6,i24gj^ ^^^^ agree (just outside errors) with previous 
data HQ- 

It has been suggested that the new measurement by Li et al. offered cleaner 



spectra with less need for subtraction o: 
this discrepancy 



background, which may possibly help to understand 
71| . Colo et al. [T^] tried to interpret the new data employing a self- 
consistent quasiparticle random-phase approximation model with Hartree-Fock-Bogoliubov 
basis, the Skyrme interaction and density dependent pairing. They reproduced the new 
GMR energies in ^^^^^^°Sn isotopes using SkM* Skyrme force and surface pairing force and 
found that the effect of pairing on Ko is very small. However, the value of Kq extracted 
from the fit to Sn isotopes is about 10% smaller than the one obtained from the fit to 
208pb, 230-240 MeV. Similar ~10% discrepancy between Kq from fits to Sn and Pb GMR 
energies has been reported in RMF calculations \(\- So, the puzzle of GMR in Sn nuclei 
remains open Keeping in mind the current unresolved situation of GMR experiments 
and theory, and that all the extracted values are likely to be rather model dependent, we 
choose the following constraint for Ko that we refer to as the SMI constraint, 

SMI : Ko = 230 ± 30 MeV (25) 

To our knowledge, there has been only one attempt to constrain the third derivative of 
the energy per particle with respect to density, the skewness K' = —Qo (see Eq. ([8])) which 
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is the next order fluctuation around the saturation density in the expansion t 

£{p)/p ^E, + {Kj2)x^ - {K'/6)x' + ■■■ . (26) 



Constraints on the skewness coefficient, which differs from author to author by a minus 



sign 



65 



77( 1 . are relatively imprecise. Farine et al. 65|] tried to find acceptable values of K' 
by self-consistent analysis of breathing mode data using a selection of Skyrme forces. They 
found a subtle correlation between Ko and K'. We adopt here their value of K' as constraint 
SM2, 

SM2 : K' = 700 ± 500 MeV. (27) 



/. [3 



The wide range of K' compensates for the fact that Farine et al. [65| used ^0=215 ±15 
MeV, 7% lower than our choice. 

Limits on the pressure-density relationship in SNM and PNM and its curvature can be 
obtained from analysis of experimental data on the motion of ejected matter in energetic 
nucleus-nucleus collisions. Recently, measurements of the particle flow in collisions of ^^^Au 
nuclei at incident kinetic energy per nucleon varyin g fro m about 0.15 to 10 GeV were anal- 



ysed The authors extrapolated available data [79| for pressure at about 2po to higher 
densities, as well as to zero temperature. The results give limits on pressure as a function 
of density (see Fig. [1]) which comprise the constraint SMS, 

SMS : P vs —{y = 0.5) =^ Flow Experiment. (28) 

Po 

Recent experimental investigation of kaon production in HIC produces a further con- 



straint on pressure at a lower density region 1.2 < p < 2.2 fm ^ 80|, |8l| (see Fig. here 
named SM4, 

SM4 : Pvs—{y = 0.5) =^ Kaons + GMR Experiments. (29) 

Po 

B. Pure Neutron Matter 

The EoS of PNM is of a particular interest, because PNM is a realistic ffist approximation 
to the baryonic matter that composes neutron stars. Most properties of neutron stars cannot 
be studied in terrestrial laboratories and theoretical models, based on effective forces, must 
be used. However, at low densities, experiments with cold Fermi atoms yield information on 
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stronKly interaction fluids, very similar to the low-density neutron matter at neutron star 



crusts 4l|]. Different density regimes can be tuned by the magnitude of the neutron Fermi 
momentum kp relative to the effective range Tq of the NN interaction in the system 
The ground state energy per particle, the EoS, can be expressed as 

-EpNM 



-^PNM 



e, (30) 



where -Epnm is the energy per particle in Eq. ([T]) with = 2^^ ^. -E'pnm = Spnm/p and 
E^^y^ = Stt^kl/IOM, with M being the mass of the nucleon. In the dilute degenerate Fermi 
gas regime, kpTo <^ 1, ^ is a constant [82]. This restricts the density below about lO^^po, 
the density at which neutrons become unbound in neutron stars. At higher densities, below 
~ O.lpo, where kpr^ ~ 1, ^ has to be replaced by a system dependent function ^(kp, ro). In 
this work we adopt the expression Ep^M /Epj^jyj by Epelbaum et al. [4J], based on next-to- 
leading order in lattice chiral effective fleld theory (NLO3), and including corrections due to 
finite scattering length, nonzero effective range, and higher order corrections, 

"f^^^ = ^ - + cikpro + C2k],m~'^ + c^klm-^ H , (31) 

PNM KpO^o 

where m,r is the pion mass. The dimensionless universal constant has been determined 
from trapped cold atom experiments with ^Li and '^''K, which yield a variety of values: 



0.32+}° y, 0.51(4) 



Jp, 0.46lJ^ [85|, and 0.39(2) |86|. Values of ^ in the literature are in 
the range 0.8 - 1.0 ( 40] and references therein). Epelbaum et al, using a simple Hamiltonian 
and only few particles in their system, took ^ = 0.31 and ^1 = 0.81 and fitted two sets of 
constants Ci, C2 and C3: set 1 (0.27, -0.44, 0.0) and set 2 (0.17, 0.0, -0.26), and obtained a 
very similar quality fits to their NLO3. We construct a constraint on energy per particle of 
PNM in the range of densities 0.01 - O.lpo shown in Fig. [3] 

c 

PNMl : -^(MeV) vs p. (32) 

Two shaded areas are based on Eq. fl3T|) with ,^1 = 0.81 and set 1 (red dashed line) or set 
2 (green dashed-dotted line). The area inside the blue solid line is based on Eq. (15CT]) . the 
unitary limit. The boundaries of all three areas are calculated by taking 0.2< ^ <0.6, which 
allows for the spread in experimental values. 

It is clear that our PNMl constraint is consistent with a relatively large range of extrap- 
olated experimental ^ values. Very recently, after our work was completed, more accurate 
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calculations and measurements were reported. The new limits on ^ are 0.37-0.38 87| and 



0.38-0.41 [88|. 1^1, related to the contact parameter in unitary Fermi gases, has been calcu- 
lated to be ^1=0.9 (jssl and references therein). These new data will be considered in future 
development of the PNMl constraint. 

New theoretical calculations also provide constrains on the EoS of low-density PNM. 
Gezerlis and Carlson ((42I and references therein) compare their EoS, obtained using the 
Quantum Monte Carlo techniques, with results of other model calculations in their Figs. 
3 and 4. Although we still keep the PNMl as selector of the low-density neutron matter 
EoS in this paper, we have constructed a band, representing the boundaries of theoretical 
predictions, in the same p/ po range used in the PNMl constraint. This band is considerably 
narrower than the band extracted from the cold atoms data. It is shown in Fig. H^, and 
merged with the PNMl constraint in Fig. |3)d. 

In the high density region, analysis of HIC data [tsI] provides a constraint on the pressure- 
density relation in PNM. 

PNM2 : P vs -^{y = 0) ^ Flow Experiment. (33) 
po 

The constraint is illustrated in Fig. [51 



C. Constraints involving both SNM and PNM 



Nuclear symmetry energy, in particular its density dependence, has received considerable 
attention during the last decade. It produces information about the isospin dependence of 
nuclear forces which is equally important in nuclear matter and finite nuclei. The current 
empirical values of J, the symmetry energy at saturation density, as predicted by different 
models, vary around the value extracted from the up-to-date finite range liquid dro ple t 

or S(p) come from several sources: HIC 



model J = 32.5 MeV 



The data 



921], Pygmy Dipole Resonances (PDR) 



91 



93 



94| . Isobaric Analog States (IAS) 



8C, 



9G- 



95| and 



numerous microscopic calculations. A systematic difference exists between predictions of 
HF (Skyrme) and RMF models 4^, l96|] spanning the range 27 < J < 38 MeV. This large 
uncertainty stems, in part, from limited experimental knowledge and understanding of the 
isospin dependence of the nucleon-nucleon interaction and, in particular, the PNM EoS. 
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Taking into account some more recent data 97|] we define this constraint as on J, 

MIXl : 30 < J < 35 MeV. 



(34) 



The density dependence of the symmetry energy, especially at super-nuclear densities, 



96 



98| and is closely relat ed to studies of 
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lOOj. In contrast to 



has direct relevance for modelling neutron stars 
neutron matter radii and the neutron skin in neutron-heavy nuclei 
the expansion of the energy per particle, in which the term containing the first derivative 
vanishes, the expansion of S{p) Eq. f lT^ contains a first order correction L at p = p^. L 
becomes an important bulk quantity that determines most of the behavior of S{p) in the 
vicinity of po- The empirical determination of L is, as for several other bulk quantities, 
indirect. The very recent constraint, based on the empirical MSL model with the MSLO 
Skyrme-like interaction, is based on analysis of isospin diffusion an d dou ble neutron/proton 



ratio in heavy-ion collisions at intermediate energies, and requires |l01 | 



MIX2 



L = 58 ± 18 MeV. 



(35) 



We comment that this value is lower than the previously accepted L = 88 ± 25 MeV 102 1 
derived by the same group but considering only the isospin diffusion data and standard 
Skjnrme forces. This lower value of L is supported by findings of Newton and Li, jl03 |. who 
used the correlation between the gravitational binding energy of a low mass neutron star 
PSR J0737-3039B and the slope of the nuclear symmetry energy at 1-2 times the nuclear 
saturation density. This correlation leads to an upper limit L < 70 M eV. It is also consistent 



with the value L = 49.9 MeV determined from the droplet model |104l | and closer to the 
most recent value from Finite Range Droplet Model (FRDM) L = 70±15 MeV 89(]. Vidana 
et al. 571] found the value of L = 66.5 MeV in their Bruckner-Hartree-Fock calculation with 
the Argonne V18 potential supplemented by a three-body force of Urbana type. As J and 
L are correlated, some investigations produce range of acceptable values for both of these 



observables (see e.g. 93|, |97|] for a recent summary). 



The next constraint involves the isospin incompressibility Kj. in Eq. f l23p for the com- 
pression modulus of a finite nucleus with mass A. This constraint on K^, which is only 
dependent on t he validity of the expansion (!36ll . from experiment was provided by Li et 
1051 1 who measured GMR strength distributions for ^^^"^^"^Sn and io6,iio-ii6qj 



al. 



70 



71 



isotopes in inelastic scattering of alpha particles. They used a simplified expression for the 
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compression modulus Ka (as compared to Eq. ([23 



= + ir,urfv4-^/=^ + Krl3^ + ^coui^, (36) 

omitting the higher order terms in A and /3. Li et al. further assumed that Ksuri = —K^oi 
and i^voi = Ko (the scahng approximation) . Kc ,m^] was taken to be (-5.2 ± 0.7) MeV. This 



value has been derived by Sagawa et al. 



106| . who investigated a microscopic structure 



of -ft^coui and its correlation with Ko, using a set of thirteen Skyrme parameterizations and 
seven RMF Lagrangians. They found this correlation rather weak and arrived to an average 
value of -ft'coui, used by Li et al. This procedure yielded the value of Kr = (-550 ± 100) 
MeV. 

It is important to realize that the K^, extracted by Li et al., includes both the volume 
and the surface components of the isospin compressibility K^-, Eq. (!24|) . It follows that, in 
order to compare microscopic model calculation with the constraint, the contribution of the 
surface-symmetry term, which is difficult to calculate e xact ly, must be at least estimated as 



well as possible and subtracted from Kr. Stone et al. |l07l | re-analysed a combined Sn+Cd 
data by Li et al. and Garg under the same conditions and found K^- = (-595 ± 154) MeV. 
Estimation of Kt-,, and K^-s from the currently available GMR data is not easy. The data 
show systematic differences, mainly dependent on methods, used for analysis by different 
groups. It is therefore necessary to include some additional constraints on the fit to obtain 
limits on Kr^v and Kr^s- Stone et al. used two assumptions. First, they required that 
Eq. ( HM holds and looked for all combinations of Kr,v and K^^s which would satisfy it. Kt^v 
is expected to be negative in line with microscopic calculations. It was varied in the region 
of -1200 < Kr,^ < MeV with K^^s in the range of -1600 < K^^s < 1600 MeV. The second 
assumption was that the expansion in terms of A~^^^ and converges at a reasonable 
rate, i.e., no higher order term are significant. They took 

< 0.5 (37) 

Simultaneous application of Eq. ( |24|) and Eq. ( 137|) . together with the assumption that K^-,^ 
is negative, allows limits on Kr^v to be extracted and taken as a constraint, 

MIX3 : -760 < K^^^ < -372 MeV. (38) 

Corresponding limits on the surface contribution to isospin incompressibility are -1110 
< Kr,s ^ 960 MeV. We comment that the condition on the ratio of the contribution of 
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the surface and volume isospin incompressibility is rather conservative. Treiner et al. 62 1 
estimated the surface contribution to the isospin incompressibihty to be ~20% of the volume 
term in ^°*Pb in the scaling approximation, implying a rapid convergence of expansion, 
Eq. fl25]) . We note that other suggested limits on K^-^^ = -370 ±120 MeV (in notation of 



the original paper [108] -ft'^'f) exist in the literatur e, bu t they are calculated, not directly 
extracted from experimental data. Patra et al. 109| estimated Kr^v and K^^s using a 
semiclassical relativistic mean- field method with interaction NLl, NL3 and NLSH. They 
obtained K^,^ = -676, -690 and -794 MeV and K^,^ = 1951, 1754 and 1716 MeV for the 
three interactions, NLl, NL3 and NLSH, respectively. 

As the next constraint, we used the known effect that at the nuclear surface, variation of 
;he d ifference between proton and neutron densities (neutron skin) is expected. Danielewicz 



110| considered this question and proposed limits to the reduction of the symmetry energy 



at Po/2 in terms of S{po/2)/ J. These limits, based on a semi-classical Thomas-Fermi model 
and a comparison of the neutron skin of ^°^Pb as calculated in this model with values from 
full mean-field models, lead to 

MIX4 : 0.57 < ^^^^ < 0.83. (39) 



Finally, Piekarewicz llll Ill2| used a parabolic approximation to the EoS and derived 



an expression for the pressure in pure neutron matter, related to the slope of the symmetry 
energy L at saturation density. Thus the "symmetry pressure" L, a quantity that influences 
the neutron-skin thickness in heavy nuclei, is directly proportional to the pressure of pure 
neutron matter as 3Ppnm{Po) / {Lp^) ~ 1 [112]. Considering the uncertainty in the number 
of terms included in the expansion, we introduce a range to this constraint as 

MIX5 :5^^^^^ = 1±0.1. (40) 

Lpo 

In terms of Skyrme model parameters, this constraint can be expressed analytically as 



'iPpNM{Po) 



yE^^^ (41) 
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where, 



T3 



_ A _^ /37r^ 
5 J 3M 
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(5 + 37) 



3^ 
5 ' 

"9v^ 



1 — 0:4) + X4 



pI" 
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g (I + X5) - -(5x5 + 4) 



hPo 



+7 



(42) 



D. Results of application of the macroscopic constraints 



Predictions of bulk properties of nuclear matter by all Skyrme parameter sets are sum- 
marized in Table HIl Table IIIII details whether a particular parameter set is consistent (+) 
or not (-) with each of the eleven constraints. 

Numerical evaluation of the compliance for individual parameterizations with the con- 
straints is given in Table IIVI The criterion for consistency with a constraint is different 
for numerical (SMI, SM2 and MIXl-5) and graphical (SM3, SM4, PNMl and PNM2) con- 
straints. For the numerical ones, we define a deviation Dev 



Dev 



Qmod Qconst 

A ' 



(43) 



where Qmod is the value of a specific quantity calculated in the model, and Qconst is the central 
value of the related constraint. The error associated with Qconst is A. If \Dev\ < 1, the model 
is consistent. In particular, for the MIXl, MIX3 and MIX4 constraints that are defined by 
a range in the form Xi < X < X2, we define the central value as Qconst = (2^2 + xi)/2, and 
the error as A = X2 - Qconst = Qconst - a^i- 

For the graphic constraints SM3, SM4, PNMl and PNM2, we consider as consistent a 
model that is inside of the corresponding band in 95% or more of the density region. In this 
case, the Skyrme interaction is rendered valid. 
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Table |V] shows the number of parameter sets (out of 240) which satisfied each con- 
straint. This table shows that no individual constraint is particularly discriminative. It is 
more interesting to consider further the parameter sets which satisfy all constraints. They 
are surprisingly few in number - only six sets LNS, NRAPR, Ska25s20, SQMC650, 
SQMC700, and SV-sym32 are selected. 

Further, considering the need to choose ranges for some constraints, we looked for sets 
that fell outside the chosen range for only one of the eleven constraints and by less then 
5% of the closest limit. This procedure yielded ten more sets GSkI, GSkll, KDEOvl, 
MSLO, Ska35s20, SKRA, SkTl, SkT2, SkT3, and Skxs20, making a total of six- 
teen Consistent Skyrme Parameterizations (hereafter CSkP) which satisfy all constraints 
on properties of nuclear matter. (We note that sets SkTla, SkT2a and SkTS a (se e Table HTj) 
have the same parameters relating to nuclear matter as SkTl, SkT2, SkT3 jll3|, and differ 
only by a choice of the spin-orbit functional, the Coulomb exchange term and the fitted 
pairing strength. SkTla, SkT2a and SkT3a are therefore not included separately in the 
analysis) . 

The values of all relevant parameters of the CSkP are given in Table I VII Table IVIII lists 
all the numerical properties of nuclear matter at saturation as calculated by the CSkP. The 
compliance of the CSkP with graphical constraints SM3, SM4, PNMl, PNM2 is illustrated 
in Figs. [M 

The range of calculated values from all the CSkP is compared with the range of each 
constraint in Table [H where, it is interesting to note, they often fall within a band much 
narrower than the imposed constraint. Figs. [TH5] show generally the same behaviour for 
constraints defined by a function. Two exceptions concern constraint PNMl (Fig. |3]) where 
a narrow band is predicted by all models except KDEOvl (high) and MSLO (low), and MIX5 
(Table IVII|) . for which all results are in the upper half of the range, although not closely 
clustered. 

As shown in Fig. Hb, where the new band from Fig. was added a to the bands and 
curves already presented in Fig. [3l two models, KDEOvl and MSLO do not satisfy the more 
stringent theoretical constraint. However, as we consider PNMl a valid constraint in the 
context of this paper, we keep the KDEOvl and MSLO parameterizations in the CSkP list. 
This issue will be revisited in future. 

The relationship between L, So=J, and pressure in pure neutron matter Pq has been 



23 



examined by Tsang et al. [90]. Based results of mass measurements, HIC, PDR and IAS, 

they produced a composite constraint on these variables which is an extended variation of 

I j 

our constraint MIX5. Fig. El taken from (^] and references therein), is the latest version 



of their Fig. 3 in 



90| . Predictions of the CSkP of this relationship all fall within the blue 



dashed rectangle, overlapping with all the constraint but showing no compatibility with the 
IAS analysis. 



IV. MICROSCOPIC AND OBSERVATIONAL CONSTRAINTS 

In addition to the eleven macroscopic constraints considered in previous sections, we 
introduce some additional, more microscopic constraints and constraints based on obser- 
vation of neutron stars and apply them only to already chosen CSkP. We find that these 
constraints significantly reduce further the number of the CSkP, eliminating GSkI, GSkll, 
MSLO, Ska25s20, Ska35s20, the SkT group, Skxs20, SQMC650, and SV-sym32 from the 
CSkP list, as discussed in detail below. 



A. The effective mass 



In non-relativistic models of the motion of a nucleon with mass M in homogeneous nuclear 
matter, the nuclear potential V(k) is momentum dependent. The concept of the effective 



mass M*, originally developed by Brueckner 114| |. leads to an equivalent description of the 
motion in which the nuclear potential V(k=0) is constant but the nucleon mass has been 
modified. It has been established that the M* is lower than the M for all potentials for 
which the low k expansion V(k)=V(0)+bk^+- ■ ■ (where b is a constant) is valid. 

This simple formalism can be extended also to momentum and density dependent poten- 
tials, such as the Skyrr ne p otential. The nucleon isoscalar effective mass M* at saturation 
density can be defined 1151 ] as 



2M? 



dE 



(44) 



where E is the energy per particle, Eq. ([T]), which leads to 45] 

M* M 

-: = ^ = (i + ^poa)-\ 



(45) 
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in terms of the Skyrme parameters. 

The isovector effective mass M* is given as 

1 



M* 

Id 



M ^ 
1 + i^PoO. 



1 + K 

where n is the enhancement factor of the Thomas-Reiche-Kuhn sum rule 
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2M 



K 



-po 



d 



d 



h? d{Tn - Tp) d{pn - Pp) 



E 



(46) 



(47) 



(notice the typographic error in 



Po 



115| Eq. (14) |ll7l |). In the above equations T,Tp and r„ are 



the total, proton and neutron kinetic energy densities and Gg and By are defined in Eq. (I^Tj) . 

In ANM, consisting of unequal amount of neutrons and protons, the nucleon effective 
mass can be written in terms of Gg and 0v as 



113 
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M. 



1_ 

M 



M 



(48) 



where f3 is the asymmetry parameter (p„ - Pp)/p and g = 1 (-1) for neutrons (protons). 
For non-standard parameterizations the Eq. (148|) becomes 

-1 



m„ 



1 + ^pQ's - (26. - Gs - U{1 + 2x,)p^ + + 2x5)p°) {3p 



(49) 



For such a case, m* and m* are obtained by making Gg — )• 6^ and 6. — )• 0. in Eqs. (H5|) 
and fj46|) . where 9^ and 9^ are defined in Eq. fl63|) . 

Constraints on the m* and m* at saturation dens ity can b e derived from experimental 
peak frequencies of giant resonances in finite nuclei ( llSl Il20l | and refere nces therein) . m* 
is solely related to giant quadrupole resonance (GQR). Kliipfel et al. |ll5| deduced an 
optim um value of m*=0.9 from the GQR in ^"^^Pb, close to the estimate of Bohias et al. 



121| m*=0.8. m* is constrained from the giant dipole resonance (GDR) which is sensitive 



to two nuclear matter variables, the symmetry energy and the enhancement factor k. Kliipfel 
et al., taking k=0.4, obtained m*=0.7, the same as the value extracted in 1221 ]. However, 
these v alues , based on the GDR in ^"^Pb, are not quite consistent with data of the GDR 
in ^^O 1231 ]. The current conclusion is that experimental GDR data on light and heavy 
nuclei cannot be satisfied simultaneously with the present form of the static HF functional. 
Additional work is needed in both theory and experiment. 

The effective mass scales the level densi ty q of single-particle (s.p.) levels in the vicinity 
of the Fermi surface g{eF) — )■ ;^g(ei?) (see 124| and references therein). This simple scaling, 
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valid in infinite matter, holds for level density of deep s.p. states in finite nuclei, but breaks 
down for s.p. states close to the Fermi level. The origin of the change in level density is 
usually seen in c oupl ing between HF s.p. modes and surface-vibration (beyond HF) RPA 
modes (see e.g. 125| ). If the coupling is taken into account, the nuclear matter effective 
mass scaling can be recovered. Such calculation is however rather complicated and the 
simplest way to fit experimental s.p. level density is to take M* being state dependent and 



equal to M at the Fermi surface 



119| . Such an approach is also necessary in fitting atomic 



mass data with conventional Skyrme forces, where a high precision fit of masses of open 
shell nuclei is not possible without a correct spacing of s.p. states close to the Fermi surface. 
However, the choice of M*=M is inconsistent, for example, in the context of formation of 
nuclear matter from isolated nuclei in neutron stars or core- collap se supernova matter. One 



pos sible remedy is to construct an extended Skyrme force [ll9l|. More recently, Satula et 



al. 



124| studied the problem of the effective mass scaling within the Skyrme energy density 



functional (EDF) method. They concluded that more detailed modelling of the two-body 
spin-orbit and tensor interaction strength reinstates the conventional m* = M* /M scaling 
and removes the inconsistency in the effective mass scaling of s.p. level densities in nuclear 
matter and finite nuclei, caused mainly by fitting strategies of the Skyrme parameters to 
incomplete Skyrme functional. 

Based on the theoretical concept of the effective mass and the experimental data on 
GQR and GDR we find a strong enough reason to eliminate all CSkP parameterizations 
with M*=M. This choice is not valid in nuclear matter (the prime concern of this work) 
and has only an auxiliary character in finite nuclei which is likely to be improved upon. The 
constraint on m* alone would eliminate the Ska and SkT family forces from the list of CSkP, 
thereby reducing their number to eleven. For the remaining, m* at saturation density is 
calculated in the range 0.603 - 0.930. Due to the weak nature of this constraint, deduced 
from experimental GDR, we do not feel there is a strong enough reason to do any further 
elimination. 

Next, we examine the density dependence of the effective neutron (Fig. [71 left panel) 
and proton (Fig. [TJ right panel) mass, m* and m*, in BEM. These are very important in 
modelling cold neutron stars. It can be seen that the set SV-sym32 yields the m* close to 
one at very low density, raising rapidly with increasing density up to about 1.7 at 3po- Such 
behaviour is not physical and is a reason for elimination of SV-sym32 from CSkP. We note 
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that some other members of the SV family 115j SV-sym28, -sym34, -K218, -K226, -K241, 
has, -kap60 and -maslO also exhibit the same feature. They all pass through a singularity 
at densities ~5po and higher. 



B. Landau parameters 



As an alternative to the Hartree-Fock approach to the properties of nuclear matter, the 
formalism of Landau theory of a normal liquid has been used. In this approach, the bulk 
properties of nuclear matter are written in terms of a two-body interaction expressed as a 
functional second derivative of the energ y per part icle with respect to occupational numbers 



at the Fermi surface. This has the form 
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= Sir. - r,)N,-' + F'^{T,.r,) + G^(a,.a,) + G',{r.,.T,){a,.a,)]PL{cos6 



(50) 



The number of states per unit energy per unit volume at the Fermi surface in symmetric 
matter is A^o = ^fs"^ where kp is the Fermi momentum. In pure neutron matter this 
quantity halves. 6 is the angle between the momenta of the interacting particles (holes). 
The sum is over angular momentum L; for the Skyrme interaction, L=0 and 1 as it contains 
only S and P wave contributions. The dimensionless parameters F and F' are directly related 
to quantities describing nuclear matter such as effective mass , incompressibility, symmetry 
energy, and the speed of sound through relationships 
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(51) 
(52) 

(53) 

(54) 



3M 1 + 1 /3Fi ■ 

The Landau parameters G and G' determine, to leading order, properties of nuclear matter 
in the spin and spin-isospin channels. We note that only six, out of the eight La ndau 



parameters in SNM, are independent, because of two Pauli principle sum rules |129| and 
conventionally F[ and G'^ are expressed as a function of the other six. In PNM, with no 
isospin degrees of freedom, only 4 parameters, Fq, Gq, Fi and Gi are non-zero. 
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It can be established that stabihty demands each of Fl, F'i^,Gl,G'i^, to be greater than 
-(2L+1) 127|, i.e., L=0 parameters must be greater than -1 and L=l terms greater than -3 
(see Figs. Eland [HI) . The most obvious justification for these conditions is the requirement that 
incompressibihty, symmetry energy (for stable HF solution for symmetric nuclear matter), 
and speed of sound be positive (see Eqs. ( 1521) . (15^ . and (IMl) ). 

Exact relation between par ameters of the Skyrme interaction and Landau parameters can 



be derived (see, e.g. |5J, Il28| ). Beside the expressions (15T]) - fl54l) for Fq, Fi, and Fq, we 



also present the remaining Landau parameters, G and G' given by 
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for SNM [52], and 



<-J0 — 'J^O "T '-^0 



/SNM 



and 



G, = Gf ^ + ^ (tl - t2 + t4/ - t5P^) 



(58) 
(59) 



for PNM 



130[ |. also valid for the non-standard Skyrme models. 



Fig. m (for SNM) and Fig. [91 (for PNM) show the results of the variation of the Landau 
parameters with density for the CSkP sets, after application of the constraint related to 
effective mass. It can be seen, at densities below about ~0.1 fm~'^, that all CSkP sets 
predict negative incompressibihty in SNM. This feature, referred to as a spinodal instability, 
should be seen as a shortcoming but also a realistic consequence of the strong correlations 
between nucleons which, at low density, cause them to form a kind of ga s-liquid separation. 
This instability was observed experimentally in heavy-ion collisions 13 1| at critical densities 
about half to two-thirds of nuclear saturation density. No such transition is predicted in PNM 
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(Fig. |9]). This instability in symmetric and BEM matt er has been consistently predicted by 



132 |. 



a variety of non-relativistic HF and RMF models 

The density dependence of the parameter Gq can be used as an indicator of a breakdown 
of spin symmetry, i.e., a transition to a spin-ordered (ferromagnetic) phase in SNM as well 
as in PNM. Such a transition could have important consequences for the evolution of th e 
proton-neutron star in core-collapse supernova and neutrino transport inside the star 133| . 
but is not as yet constrained experimentally. Theoretical studies of the spin-ordered phase 
yield rather contradictory results. Skyrme int eractions predict such a transition at low 



densities (below 3po) in PNM, SNM and ANM 133|-137|. Relativistic DBHF calculations 
with an effe ctive Lagrangian also predict a transition to a spin-ordered phase at several times 



Po (see 



138l | and references therein). On the contr ary, realis tic NN interactions sup press such 



a transition up to high densities in BHF models 



lowest order constrained variational method 
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140j, the AFDMC method [l4l|, and 



As shown in Fig. [HI the transition to spin-ordered matter is predicted in SNM at densities 
below Po for the GSkI and GSkll, and at around 1.5 po for MSLO, Skxs20 and SQMC650 
parameterizations. In PNM, five parameterizations predict the transition at densities below 
and around p^, GSkI, GSkll, MSLO, Skxs20 and SQMC650 and two at around 1.5 po, SKRA 
and SQMC700. 

These features are not realistic and we eliminate GSkI, GSkll, MSLO, Skxs20 and 
SQMC650 in their present form from the CSkP list as they did not satisfy both, the SNM and 

56|, the spin-spin 



PNM constraint. However, as demonstrated by Margueron and Sagawa 
and spin-isospin instabilities can be removed if additional density dependent terms, affecting 
only the spin and spin-isospin channels, are included in the standard Skyrme Hamiltonian. 
The contribution of the new terms to the mean field is zero for spin-saturated systems. 
Consequently the properties of the original Skyrme interaction are not changed in this case. 
However, in nuclear matter four new parameters have to be adjusted to values of Landau 
parameters at saturation density extracted from the G-matrix. This procedure is rather 
involved and is beyond the scope of the present work. We suggest that parameterizations 



56| 



m 



SKRA and SQMC700 would be the best candidates for application of treatment 
future. 

We note that the density dependence of the parameter Gq, indicates a spin-isospin insta- 
bility if Gq falls below -1. Such instability has been interpreted as the appearance of a pion 
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condensate 1431 . Il44l |. Of the remaining parameterizations, only KDEOvl and LNS predicts 
such transition in SNM below 3po- AH the other Landau parameters are within the natural 
constraints. 



C. Density dependence of the symmetry energy 

One rather surprising result, which came out of our analysis, is that the CSkP exhibit 
a growing spread in density dependence of the symmetry energy beyond about twice the 
nuclear saturation density. This feature is illustrated in Fig. [TU] and Table rvTTTl and suggests 
that constraining the derivatives of the symmetry energy at the saturation point is not 
sufficient for controlling the slope of S{p) at higher densities. Clearly, more experimental 
data is needed to constrain the Skyrme interaction at super-saturation densities. 

It turns out that, considering the symmetry energy being the difference between the 
energy per particle in pure neutron and symmetric matter (to the first order), the factor 
which mainly determines the behavior of the symmetry energy with increasing density is 
the pure neutron matter EoS. In Fig. [11] we see energy per particle in PNM and SNM as a 
function of density as calculated with Skxs20, QMC700 and GSkll parameterizations. These 
forces were chosen as they correspond to the top, middle and bottom curves in Fig. [TOl We 
see clearly that the energy per particle for SNM are rather similar in all three panels, but 
for PNM it exhibits systematically a different pattern. In a sense this is not surprising. 
Skyrme parameterizations are usually fitted to properties of nuclei with either N = Z or 
a low value of isospin. The EoS for PNM is well constrained at low densities; at super- 
saturation densities we have to rely on theoretical models or seek indirect evidence from 
astrophysical extrapolations, e.g. to neutron stars. So constraining the PNM EoS by study 
of very neutron rich heavy nuclei should be desirable. 

D. High-mass cold neutron stars 

One possibility is to use the Skyrme EoS in cold neutron star models up to 3po- Here 
the Skyrme interaction is applied to n + p + e + fi BEM rather then symmetric or pure 
neutron matter. There are no constraints available from terrestrial experiments at present, 
as the heavy-ion reactions are too fast to build equilibrium conditions with respect to weak 
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interactions such asp + e ^ p + jj, . 

The findings of Ref. f98j indicate that, if one accepts vahdity of the Skyrme interaction 



at densities up to 10 times nuclear saturation density, only parameterizations predicting 
growing (or monotonously slowly decreasing - for details see Ref. josi]) symmetry energy 
with increasing density can be used to generate stable neutron sta r models with mass and 
radius consistent with currently available observational data 145| . We observe the same 
phenomenon here for selected CSkP as illustrated in the left panel of Fig. [121 Here the 
mass-radius relation for neutron star models was calculated using a BEM Skyrme EoS; 
since only a part of the neu tron star is in the BEM phase, we used the Baym-Pethick- 
Sutherland EoS (BPS) 146l | at lower densities matching the Skyrme based part at p ~ 
0.1 fm~^ and going down to p ~ 6.0x10^^^ fm~'^. The observed spread in the maximum 
mass models is not unexpected; it is related to the different extrapolation properties of 
the Skyrme interaction to densities well beyond the vali dity of the Skyrme model. A very 
recent observational finding of a massive neutron star 1471 ] with Mg=1.97±0.04 M© and 
the central density less than ^ lOpo provides a strong constraint on EoS of BEM. This 
constraint would certainly eliminate SKRA, SQMC700, and LNS EoS which predict lower 
maximum mass models. Moreover, the central densities of all maximum mass neutron star 
models predicted by CSkP, including those within the window set by the results of Demorest 
et al. for gravitational mass, predict higher central densities than is allowed (in the region 
for ll-13po). We thus conclude that extrapolation of the Skyrme model beyond its validity 
range to high densities does not predict cold neutron stars in agreement with the recent 
observation. We do not eliminate any Skyrme interaction from the CSkP list on the basis of 
this constraint because it requires extension of the Skyrme model outside its validity range. 

An alternative is to use the Skyrme interaction within its validity range argued in this 
paper (up to about 3po) to construct the EoS of BEM and match it to an established high 
density EoS, as well as with B PS E oS at lower densities. Such an EoS was usually thought 



to be the Bethe- Johnson EoS 



148| based of the Reid potential and including hyperons in 



the composite matter at high densities. However, the Demorest et al. observation rules 
out this EoS as it predicts maximum mass of the neutron star models to be only ^ 1.85 
Mq with centr al density ~ 12po. Therefore we use the Full-Quark-Meson-Coupling model 
(FQMC) 149l | which includes full baryon octet in the high density matter and provides 
high mass neutron star models in agreement with observation, as shown in the right panel 
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of Fig. [T2I The maximum mass is clearly determined by the FQMC model and only the 
variations in radii of smaller neutron star models with smaller central densities are due to 
different Skyrme EoS. Observational data on neutron star radii are very poor at present. 
However, new observational techniques are being developed and radii may be known within 
a few percent in the near future. They could be then used as a useful constraint on the 
Skyrme parameterizations performance in a high density neutron rich environment. 



E. Low-mass cold neutron star 



Observation of the double pulsar J0737-3039 provides another, s tring ent, constraint on 
the neutron star EoS and its interpretation by Podsiadlowski et al. 150[, and hence on the 
effective nucleon-nucleon interaction in stellar matter. This constraint is important, as the 
central density of pulsar B is only 2-3po and thus the Skyrme interaction is expected to be 
fully applicable. 

The constraint concerns the ratio of the gravitational mass of pulsar B to its baryonic 
mass. The gravitational mass of Pulsar B is very precisely known Mg=1.249±0.001 M©. 
Estimates of the baryonic mass depend upon its detailed mode of formation. If the pulsar 
B was formed from a white dwarf with an 0-Ne-Mg core in an electron-capture supernova, 
assuming no or negligible loss of baryonic mass during the collapse, the newly born neutron 
star will have the same baryonic mass as the precollapse core of the progenitor star. As 
modelled by Podsiadlowski et al. 150l |. the baryonic mass of pulsar B is then between 
1.366 - 1.375 Mq. The range reflects uncertainties in modelling of the progenitor core such 
as electron-capture rates, nuclear network calculations and Coulomb and general relativity 



corrections. Another simulation of the same process by Kitaura et al. |151| gave Mq 
1.360±0.002 Mq where the range of the result was mainly due to uncertainty in the EoS 
and the estimated small mass loss during the collapse. 

For any neutron star matter EoS the relation between the gravitational and baryonic 
mass can be calculated. Fig. [13] shows the results for the remaining five CSkP to be checked 
against narrow windows given by the two models of pulsar B. All of the CSkP but NRAPR, 
which is slightly shifted, agree remarkably well with the result of Kitaura et al. 15l| thus 
supporting the concept of some baryonic mass loss during the collapse. 



32 



V. DISCUSSION AND SUMMARY 



In embarking on this project we hoped the present work might lead to better under- 
standing of how a universal Skyrme force might best be achieved. In this section we discuss 
the degree to which this hope has been accomplished. A distinction was made between 
standard and non-standard forms of the Skyrme interaction, the latter including additional 
terms added relatively recently. The overall results are discussed below for the two forms. 



A. Standard parameterizations 

Only six of the over 224 standard sets satisfied all the constraints. Examination of 
Table IIIII reveals that some families of sets show systematic patterns of failure, as discussed 
briefly below. Other parameterizations, among them several commonly used for years in 
modelling finite nuclei, for example the SIII, SkM*, SkP, SGI, SGII and the Ski and SKX 
families, each fail different constraints with no apparent pattern. 

This very variable performance makes it rather difficult to identify a particular term or 
terms in the Skyrme energy functional as responsible. Among the more systematic inconsis- 
tencies, BSkl-17 forces, with the exception of BSkl4, do not satisfy constraints PNM2, MIX2 
and MIX3, a feature shared with most of their predecessors, the SkSC and MSk families. 
They predict too low pressure in PNM as a function of increasing super-saturation density. 
They underpredict both the symmetry energy and its derivative at the saturation density 
(see Table HTl and Fig. [6]). In addition, the failure of MIX3 indicates too high volume part 
of the isospin incompressibility, i.e. they overpredict isospin dependence of the curvature of 
the EoS of ANM E(p,/3) at the saturation density p(/3) (see Eq.([I7])). 

None of the well known Lyon forces satisfy constraint MIX3, although they all pass 
constraint MIX2 (except for SLylO), similarly to the BSk forces. It is interesting to note 
that SLy230a force [45;], especially developed for modelling of neutron stars, fails constraint 
PNMl, i.e. does not have the correct density dependence of pressure at sub-nuclear densities, 
especially important in modelling of neutron star crusts. The ov er-pr ediction of Kt-^v is also a 



152| using a fitting protocol 



feature of recently developed extensive set of new Skyrme forces 
similar to that used for construction of Lyon forces. 

Failure of MIX3 may seem to be a minor defect but it is persistent and points to the 
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isospin part of the Skyrme force. Isospin dependence of the curvature of the EoS of ANM 



plays an important role in modelling of giant resonances and heavy -ion collisions |96l 



108 | 



The very recently developed modified Skyrme-Like (MSL) model [lQl|, which is expressed 
in terms of 9 macroscopic observables that are either constrained experimentally or well 
known empirically, offers another non-traditional approach to the construction of a Skyrme 
parameterization. It expresses the standard Skyrme parameters in terms of these observables 
and provides a parameterization MSLO that complies with all but one of the constraints 
studied in this work. It predicts spin instability around nuclear saturation density, which 
would be a problem, especially in using this force to model neutron stars. 

The successful CSkP sets, several of which were unfortunately infrequently used, do not 
share much common ground. Their individual parameters, listed in Table IVTI are too spread 
to give a useful guidance for construction of a more general "consistent" set. 

In recent years some effort, aimed at giving more physical relevance to particular terms in 
the Skyrme energy functional, has lead to its re-expression in terms of the coupling constants 
(some of them density dependent) involving linear combinations of the individual parameters 
in Eq. ([T]). Four of these coupling constants, relevant for calculations of nuclear matter, are 

3 3 

CP = -to + —hp':, 

8 48 

c: = —ti + -h(- + x2) 

16 4 4 

CI = -lh{^ + x,) + h,{^ + x,). (60) 

Numerical values of these coupling constants at saturation density are given for CSkP in 
Table IIXI Again, unfortunately, they do not exhibit any apparent regularity. 



Another combination of parameters, involving only ti,t2,xi,X2, has been introduced 45 1 

Bs = 3ti + (5 + 4x2)^2, 

= ti(2 + si) +t2(2 + a;2), 

6sym = 3tia;i - t2(4 + 5a;2), 

= h{l-Xi) + 3t2il + X2). (61) 

6s and 6v were used in Sec. IIV Al in connection with the effective mass. Ogym is used in 
calculation of the symmetry energy and Bn appears in the expression for the EoS of pure 
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neutron matter. Numerical values, given in Table |IX] for the CSkP, again show a large 
scatter. 

We are forced to the conclusion that our analysis of the performance of the standard 
parameterizations in nuclear matter does not offer any clear direction for the development 
of a unified, generally applicable, Skyrme parameterization. Indeed, the overwhelming im- 
pression, however well intentioned, is that they are merely empirical attempts to describe 
nuclear matter related phenomena. We see the problem as lying first with the lack of a 
direct connection between the terms of the Skyrme energy functional and specific physics 
observables, and second with the strong correlations between the parameters. 



B. Non-standard parameterizations 

In addition to the standard definition of the Skyrme energy density functional, some 
recently developed Skyrme models include higher order terms, thus introducing additional 
parameters. The standard form (for application in nuclear matter) depends on 9 parameters, 
tj, Xj (with j = — 3) and cr, whereas the extended Skyrme models considered in this work 
have t^i and crsj (with i = 2,3), and in addition, totaling 15 adjustable parameters. 

Out of all Skyrme models we analysed, 16 are non-standard, and some of those share 
their origin with the standard ones and are closely related to them. For example, the BSk 
family consists of 22 individual models, 4 of them (BSkl8-21) being non-standard. The 
BSkl8 model behaves as a standard at SNM, since the contribution of the additional terms 
included adds up to zero. The extended BSk parameterizations were generated by Goriely 
et al. to improve calculation of nuclear masses at the neutron drip-line for Z, > 8 and 
Z < 120. The terms containin g ^4 and were introduced to ensure a better description for 



homogeneous neutron matter 



1301 1 . With the same aim, and by adding to the conventional 



Skyrme forces higher order density term s in the EoS (^32 and contributions of Eq. ([T])), 



Farine et al. suggested the SkPS.l force 



119| . These authors claim that this force fitts well 



the nuclear masses of spherical nuclei, and may be useful to describe stelar collapse processes. 
Following the same protocol of the "Saclay-Lyon" group (SLy-forces), and by using a better 
control regarding the spin-isospin instabilities via Landau parameters, T. Lesinski et al. 



153| have also developed three non-standard forces (f_, fo and f+). They have attempted to 



constrain the effective neutron mass in such a way that m* < m* m* = m* and m* > m* 
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respectively for f_, fo, and f+. By keeping the term in the EoS, Eq. ([T]), and aiming 
at a good description of excited states of finite nucl ei, K. Kr ewald et al. suggested six 



new non-standard Skyrme parameterizations: GSl-6 [15J, Il55l ]. Agrawal et. al. [51], by 
exploring the extended density-dependent Skyrme effective forces for normal and isospin-rich 
nuclei for neutron stars, parameterized two non-standard Skyrme forces (GSkI and GSkll) 
by adjusting t^2 and ^33. These models were able to fit consistently thirteen finite nuclei:^^0, 
2^0, i^Ca, ^^Ca, ^^Ni, ^SNi, ^^Ki, ^^Ni, ^ssr, ^^Zr, ^oosn, i^ag^^ ^nd ^ospb. The breathing 
modes for ^°Zr and ^°'^Pb were also well described. 

From all the above listed Skyrme non-standard forces, only two, namely, GSkI and GSkll, 
satisfied the macroscopic constraints but failed the microscopic ones, namely the value of the 
Landau parameter Gq. Inclusion of any nonstandard piece in the energy density functional, 
Eq. (II]), inevitably affects parameters of the standard part as both contributions have to 
be compensated to fit experimental data. Nevertheless, it is instructive to investigate the 
non-standard contribution to the energy per particle of symmetric matter, obtained from 
Eq. ([T]) and shown in Fig. [HI We see that non-standard terms may either increase repulsion 
(attraction) by a positive (negative) term in the effective Skyrme force. It is interesting 
to notice that both GSkI and GSkll forces receive very similar large negative contribution 
from the non-standard terms apparently needed to compensate repulsion coming from the 
standard part of the interaction. However, because of this delicate balance between the 
standard and non-standard terms it is difficult to find any general trend. 

As in the case of the standard Skyrme parameterizations, one can also define the following 
set of couphng constants {43 1, 

3 

48 



1=1 

3 



16 4 V4 J 16' " 4 V4 
written as a function of the non-standard parameters. Notice that these equations are 
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generalized forms of the Eq. (160|) . The same occurs for the quantities 45 1 



sym 



e 




n 



es + 3t4p' + t5(5 + 4x5)p^ 

+ Ui2 + x4)p^ + ^5(2 + xr,)p\ 

©sym + ^hXip^ -h{4: + 5X5 

en + t4(l-a;4)p^ + 3t5(l+a;5)p^. 



(63) 



All these values are shown in Table IIXI 

Although it may be useful to extend the Skyrme functional to improve results in particular 
physical situation, it does not seem to be a way forward to finding a recipe for getting the 
Skyrme model under control. 

VI. CONCLUSIONS 

We have examined the performance of 240 different Skyrme model parameterizations 
in nuclear matter by comparing their predictions of behaviour in eleven areas in which 
experimentally/empirically derived constraints exist. The chosen macroscopic constraints 
cover a wide range of properties related to Symmetric Nuclear Matter (SM1-SM4), Pure 
Neutron Matter (PNM1-PNM2) and both SNM and PNM (MIX1-MIX5). Of the Skyrme 
models six satisfy all the constraints whilst 66 satisfy all but one. For ten of the 66 the single 
failure is narrow, less then 5%. Including these yield a final list of 16 consistent models, 
the CSkP set GSkI, GSkll, KDEOvl, LNS, MSLO, NRAPR, Ska25s20, Ska35s20, 
SKRA, Skxs20, SQMC650, SQMC700, SkTl, SkT2, SkT3, and SV-sym32. The 
parameters of these interactions are summarized in Table |Xl 

As an additional step, we considered four microscopic constraints arising from giant 
resonance experiments on nucleon effective mass. Landau parameters in SNM and PNM, and 
observation of low- mass neutron stars. With these microscopic constraints taken into account 
the successful set reduces to five, KDEOvl, LNS, NRAPR, SKRA, and SQMC700, 
the CSkP* set. 

A fifth microscopic constraint, maximum mass and the corresponding central density of 
high-mass neutron stars creates a fundamental obstacle to applying Skyrme (nucleon only) 
models in neutron star modelling since it requires extrapolation to densities above the range 
of validity. None of the CSkP models produce a maximum mass neutron star model with 
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central density in line with observation. Thus if this constraint were applied, all CSkP would 
fail. 

We were unable to identify any regularities, either in single parameters or in their com- 
binations, to identify a unique quality of the CSkP* sets. This is hardly surprising when 
we consider the number of parameters and their correlations, and it is tempting to suggest 
that some of the sets satisfied all the constraints by a fortunate accident. This looks partic- 
ularly likely in the case of the KDE family. Of very close KDEO parameter sets (KDEOvO, 
KDEOvl), which fit the same experimental data but differ only by starting conditions for the 
simulated-annealing fit procedure, only KDEOvl passes our constraints. It may be signifi- 
cant that in this procedure particular attention was paid to the inclusion of nuclear matter 
quantities in the fit. Also energies of the giant monopole resonances were included directly 
to the fit for the first time. 

The example of KDEOvl indicates that the inclusion as many constraints as possible, both 
macroscopic and microscopic, in the fitting protocols of the Skyrme interaction is essential. 
For example, we believe that the symmetry energy plays a key role in the behaviour of 
nuclear matter. Therefore the correct determination of the PNM EoS is imperative. Whilst 
there are no direct experimental or observational data on PNM at super-saturation densities, 
promising ab initio theoretical predictions and indirect experimental data from cold atoms 
at sub-saturation densities are becoming available. Any further development in this area is 
very desirable. 

It may be also revealing that in construction of three out of five CSkP* parameterizations, 
LNS, NRAPR, and SKRA, the EoS of nuclear matter, used in the fit, was derived from 
realistic potentials. Such an approach amplifies the role of the microscopic physics input in 
the effective Skyrme interaction model. The usual practice of including basic properties of 
nuclear matter only at saturation density is not sufficient. The density dependence of these 
observables (within the range of applicability of the Skyrme model), which may considerably 
infiuence the fits must be included. 

The last parameterization on the CSkP* l ist, SQMC700 , is unique in its derivation from 



the Quark- Meson-Coupling model (FQMC) 1491 . 1 1561 . Il57l |. which includes the full baryon 



octet at high densities and is relativistic. The structure of the QMC Hamiltonian and its 
Skyrme equivalent (FQMC limit at low densities) are somewhat different, which shows up 
most markedly in the difference between the values of the incompressibility for FQMC EoS 
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149| and QMC forces considered in this paper. Nevertheless, the fact that in FQMC the 



many-body interaction is directly related to the response of quark structure to the nuclear 
environment may cause the very simple QMC Skyrme parameterization, with a;i=a;2=a;3=0 
and cr=l/6, to be consistent with all the constraints. 

It is important to stress that our present work is restricted to examining the performance 
of the Skyrme interaction in nuclear matter related scenarios. It is generally true that only 
a limited effort was spent to make the CSkP* perform well in wide-ranging application in 
finite nuclei when they were derived in the original papers. In some sense, this may be a 
positive factor as it seems obvious that Skyrme interactions, constructed with emphasis on 
nuclear matter properties, will do better on constraints derived from nuclear matter then 
interactions heavily biased towards properties of finite nuclei. As a follow-up project of this 



work, the performance of CSkP in finite nuclei will be studied in more detail |158| . It will 
be interesting to test the quality of the Skyrme functional by extending the fitting protocol 
to include not only all the constraints studied in this work but also requirements based on 
the most up-to-date finite nuclei data, including drip-line and superheavy nuclei. It may 
turn out that the ambition to fit using such a de taile d protocol may be asking too much of 
the Skyrme model (with its known deficiencies 123l |) but the attempt may also lead to a 
parameterization with increased predictive power. It is our opinion that a parameterization, 
successful in reproducing only a selected subset of available data but failing the rest, does 
not have credible predictive power and does not progress the field. 

Results obtained in this paper should be seen as the first step in a global effort to find 
the best possible Skyrme interactions for use in modeling of nuclear matter. We intend to 
monitor developments both in keeping up-to-date the existing constraints a nd ad opting new 
ones. Progress in ab-initio calculations of inhomogeneous neutron matter jlSQ^ and chiral 



effective field theory 15|, ll6| are good examples of sources of new constraints which will be 



taken into account in future. New experimental results on giant resonances, neutron skin in 
heavy nuclei and heavy ion collision, as well as new astrophysical observations will further 
improve the set of benchmark constraints, which may shed more light on the structure and 
applicability of the Skyrme interaction. Should new Skyrme parameter sets appear in the 
literature, we intend to catalogue and test them using the most complete set of constraints 
available to us. 

The outcome of our analysis of all, standard and non-standard, parameterizations does 
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not offer a final solution to the "Skyrme proliferation" problem. Neither does it provide gen- 
eral guidance for construction of more Skyrme parameter sets. Production of new parameter 
sets having limited range of application should not be encouraged. Rather, more emphasis 
should be put on better understanding of the existing models, which should be further tested 
against an extended number of refined constraints including both, nuclear matter and finite 
nuclei related properties, with equal emphasis. The Skyrme interaction has played a domi- 
nant role in low energy nuclear physics for decades. The approach suggested can lead to the 
final judgment whether or not this interaction includes enough physics not only for a suc- 
cessful interpretation, but also for a prediction, of the rich variety of data and observations 
on nuclear and astrophysical systems available today and expected in future. 
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TABLE I: List of macroscopic constraints and the range of their experimental/empirical values, 
density region in which they are vahd and the corresponding range as found using successful Skyrme 



parameterizations (CSkP). For more explanation see text. 



Constraint 


Quantity 


Eq. 


Density Region 


Range of constraint 
exp / emp 


Range of constraint 
from CSkP 


Ref. 


SMI 


Ko 


10, 113 


po (fm-3) 


200 - 260 MeV 


202.0 - 240.3 MeV 


m 


SM2 


K' = -Qo 


©.Cell 


Po (fm-3) 


200 - 1200 MeV 


362.5 - 425.6 MeV 


[65] 


SM3 


P{P) 


® 


2 < < 3 

Pa 


Band Region 


see Fig. [T] 


[78] 


SM4 


P(P) 


m 


1.2 < < 2.2 


Band Region 


see Fig. |2] 


[80] 


PNMl 


EpNM 
^PNM 




0.014 < < 0.106 

Po 


Band Region 


see Fig. |3] 


[39, 40] 


PNM2 


P{p) 




2 < < 3 

Po 


Band Region 


see Fig. [5] 


[78] 


MIXl 


J 


m 


Po (fm-3) 


30 - 35 MeV 


30.0 - 35.5 MeV 


[44] 


MIX2 


L 


Go)) 


Po (fm-3) 


40 - 76 MeV 


48.6 - 67.1 MeV 


[lOlJ 


MIX3 






Po (fm-3) 


-760 - -372 MeV 


-407.1 - -360.1 MeV 


[107] 


MIX4 


S(pa/2) 
J 




Po (fm-3) 


0.57 - 0.86 


0.61 - 0.67 


[110] 


MIX5 


^PpNM 


m 


Po (fm-3) 


0.90 - 1.10 


1.02 - 1.10 


[112] 
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TABLE II: Saturation properties of all Skyrme parameterization used in this work. All entries are in MeV, except for the saturation density po in fm ^ 
and the dimensionless effective mass m* = M* /M. 



Model 




Po 


E 




K 




K' 


1 


J 




L 




-^sym 


Qsym 




m 




BSkl [1601 


0, 


.157 


-15, 


,81 


231, 


,31 


385, 


,59 


27, 


,81 


7, 


,19 


-281 


.83 


606. 


,46 


-312, 


,97 


1, 


,05 


BSk2 [1611 


0, 


.157 


-15, 


,80 


233, 


,65 


380, 


,07 


28, 


,00 


7, 


,98 


-296, 


.98 


557. 


,91 


-331, 


,87 


1, 


,04 


BSk2' fl61] 


0, 


.157 


-15, 


,79 


233, 


,32 


380, 


,86 


28, 


,00 


7, 


,79 


-298, 


.02 


558. 


,62 


-332, 


,04 


1, 


,05 


BSk3 [162] 


0, 


.157 


-15, 


,81 


234, 


,81 


380, 


,83 


27, 


,93 


6, 


,78 


-306 


.90 


550, 


,34 


-336, 


,59 


1, 


12 


BSk4 [163] 


0, 


.157 


-15, 


,77 


236, 


,84 


367, 


,17 


28, 


,00 


12, 


,54 


-265, 


.93 


558, 


,40 


-321, 


,74 


0, 


,92 


BSk5 [163] 


0, 


.157 


-15, 


,80 


237, 


,19 


367, 


,86 


28, 


,70 


21, 


,41 


-240 


.30 


499, 


,92 


-335, 


,56 


0, 


,92 


BSka [163] 


0, 


.157 


-15, 


,75 


229, 


.14 


370, 


,64 


28, 


,00 


16, 


,84 


-215, 


.19 


603, 


,54 


-289, 


,01 


0, 


,80 


BSk7 [163] 


0, 


.157 


-15, 


,76 


229, 


.26 


370, 


,92 


28, 


,00 


17, 


,99 


-209, 


.35 


598, 


16 


-288, 


,18 


0, 


,80 


BSk8 [164] 


0, 


.159 


-15, 


,83 


230, 


.31 


372, 


,39 


28, 


,00 


14, 


,85 


-220 


.88 


624, 


,89 


-285, 


,98 


0, 


,80 


BSk9 [165] 


0, 


.159 


-15, 


,92 


231, 


.32 


374, 


,67 


30, 


,00 


38, 


,29 


-153, 


.70 


482, 


,61 


-321, 


,44 


0, 


,80 


BSklO [166] 


0, 


.159 


-15, 


,91 


238, 


.83 


370, 


,34 


30, 


,00 


37, 


,24 


-194, 


.90 


396, 


,99 


-360, 


,60 


0, 


,92 


BSkll [166] 


0, 


.159 


-15, 


,86 


238, 


.09 


369, 


,18 


30, 


,00 


38, 


,36 


-189, 


.81 


390, 


14 


-360, 


,48 


0, 


,92 


BSkl2 fl66] 


0, 


.159 


-15, 


,86 


238, 


,06 


369, 


,11 


30, 


,00 


38, 


,01 


-191, 


.35 


392, 


,53 


-360, 


,47 


0, 


,92 


BSkl3 [166] 


0, 


.159 


-15, 


,86 


238, 


,09 


369, 


,17 


30, 


,00 


38, 


,82 


-187, 


.90 


386, 


,57 


-360, 


,65 


0, 


,92 


BSkl4 [167] 


0, 


.159 


-15, 


,85 


239, 


,33 


358, 


,67 


30, 


,00 


43, 


,91 


-152, 


.02 


388, 


,27 


-349, 


,68 


0, 


,80 


BSkl5 [168] 


0, 


.159 


-16, 


,04 


241, 


,56 


363, 


.14 


30, 


,00 


33, 


,60 


-194, 


.35 


466, 


,51 


-345, 


,43 


0, 


,80 


BSkl6 [169] 


0, 


.159 


-16, 


,05 


241, 


,67 


363, 


.58 


30, 


,00 


34, 


,88 


-187, 


.37 


461, 


,89 


-344, 


,17 


0, 


.80 


BSkl7 [170| 


0, 


.159 


-16, 


,06 


241, 


,69 


363, 


.62 


30, 


,00 


36, 


,29 


-181, 


.84 


450, 


,48 


-344, 


,97 


0, 


.80 


BSklS [52] 


0, 


.159 


-16, 


,06 


241, 


,79 


363, 


.82 


30, 


,00 


36, 


,22 


-180, 


.90 


454, 


,52 


-343, 


,71 


0, 


,80 


BSkl9 [130] 


0, 


.160 


-16, 


,08 


237, 


,33 


297, 


.89 


30, 


,00 


31, 


,90 


-191, 


.44 


472, 


,94 


-342, 


,79 


0, 


,80 


BSk20 [130| 


0, 


.160 


-16, 


,08 


241, 


,39 


282, 


,26 


30, 


,00 


37, 


,38 


-136, 


.49 


549, 


,73 


-317, 


,05 


0, 


,80 


BSk21 [130] 


0, 


.158 


-16, 


,05 


245, 


,80 


274, 


,09 


30, 


,00 


46, 


,56 


-37 


.20 


709, 


,66 


-264, 


,62 


0, 


,80 


E [171] 


0, 


.159 


-16, 


,13 


333, 


,46 


63, 


,72 


27, 


,66 


-31, 


,27 


-570 


.73 


448, 


,61 


-389, 


,09 


0, 


,87 


Es [171] 


0, 


.163 


-16, 


,02 


248, 


,60 


352, 


,41 


26, 


.44 


-36, 


,86 


-457 


.76 


880, 


,01 


-288, 


,86 


0, 


,84 


f- [153] 


0, 


.162 


-16, 


,02 


230, 


,01 


404, 


,93 


32, 


.00 


43, 


,78 


-105, 


.08 


654, 


,90 


-290, 


,70 


0, 


,70 


f+ [153] 


0, 


.162 


-16, 


,04 


230, 


,01 


406, 


,17 


32, 


.00 


41, 


,54 


-117, 


.98 


661, 


,07 


-293, 


,85 


0, 


,70 


fO [153,1 


0, 


.162 


-16, 


,03 


230, 


,01 


405, 


,45 


32, 


.00 


42, 


,41 


-113, 


.41 


657, 


,36 


-293, 


,11 


0, 


,70 


FPLyon [172] 


0, 


.162 


-15, 


,92 


217, 


,03 


399, 


,45 


30, 


.93 


42, 


,76 


-135 


.60 


485, 


,79 


-313, 


,47 


0, 


,84 


Gs [171] 


0, 


.158 


-15, 


,59 


237, 


,29 


348, 


,79 


31, 


.13 


93, 


,31 


14, 


,07 


-26, 


,92 


-408, 


,61 


0, 


.78 


GSl [154] 


0, 


.159 


-16, 


,03 


235, 


,09 


812, 


,19 


28, 


.86 


50, 


,22 


-58 


.00 


965, 


,37 


-185, 


,83 


0, 


,60 


GS2 [m] 


0, 


.159 


-16, 


,01 


300, 


,14 


321, 


,65 


25, 


.96 


30, 


,27 


-188, 


,78 


467, 


,72 


-337, 


,95 


0, 


,60 


GS3 [154] 


0, 


.159 


-16, 


,00 


399, 


,91 


-428, 


,61 


21, 


,49 


-0, 


,44 


-389, 


,11 


-293, 


,19 


-386, 


,02 


0, 


,60 




u. 


. loo 


-10, 


.yo 


ZOO, 


1 

.10 


o40. 


■ OO 




.OO 


1 Q 

-lo 


, ( U 


-161 


.42 


1121, 


,05 


-116, 


,53 


u. 




GS5 [154] 


0, 


.158 


-15, 


,91 


299, 


.20 


358, 


,53 


18, 


.70 


-12, 


,14 


-290, 


.66 


625, 


,40 


-232, 


,35 


0, 


,80 


GS6 [154] 


0, 


.159 


-16, 


.04 


400, 


.86 


-383, 


,63 


14, 


,33 


-42, 


,98 


-492, 


.21 


-125, 


,65 


-193, 


,21 


0, 


,80 


GSkI [51] 


0, 


.159 


-16, 


.02 


230, 


.21 


405, 


,58 


32, 


.03 


63, 


,45 


-95 


.29 


293, 


,44 


-364, 


,19 


0, 


,78 


GSkll [51] 


0, 


.159 


-16, 


.12 


233, 


.40 


398, 


,73 


30, 


.49 


48, 


,63 


-157 


.83 


310, 


,27 


-366, 


,54 


0, 


,79 


KDE [173] 


0, 


.164 


-15, 


.99 


223, 


.90 


381, 


,81 


31, 


.97 


41, 


,42 


-141, 


.83 


543, 


,33 


-319, 


,71 


0, 


,76 


KDEOv [173] 


0, 


.161 


-16, 


.10 


228, 


.71 


373, 


,39 


32, 


.98 


45, 


.21 


-144 


.78 


523, 


,27 


-342, 


,24 


0, 


,72 


KDEOvl [173] 


0, 


.165 


-16, 


.23 


227, 


.54 


384, 


,86 


34, 


.58 


54, 


.69 


-127 


.12 


484, 


,45 


-362, 


,78 


0, 


,74 


LNS [118.1 


0, 


.175 


-15, 


.32 


210, 


.78 


382, 


,55 


33, 


.43 


61, 


.45 


-127, 


.36 


302, 


,46 


-384, 


,55 


0, 


,83 



Continued on next page 
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TABLE II - continued from previous page 



Model 




Po 


E 





K 





K' 


1 


J 




L 




-^sym 


Qsyi 


in 






m* 


MSkl [174] 


0. 


,157 


-15, 


,83 


233. 


.73 


379. 


.97 


30. 


,00 


33. 


,92 


-200. 


,02 


448, 


.66 


-348. 


,39 


1.00 


MSk2 [1741 


0. 


,157 


-15, 


,83 


231. 


.65 


386. 


.21 


30. 


,00 


33. 


,35 


-203. 


,44 


449, 


.71 


-347. 


,94 


1.05 


MSk3 [174] 


0. 


,157 


-15. 


,79 


233. 


,25 


379. 


,01 


28. 


,00 


7, 


.04 


-283. 


,52 


615 


.65 


-314. 


,33 


1.00 


MSk4 [174] 


0. 


,157 


-15, 


,79 


231. 


.17 


385. 


.26 


28. 


,00 


7, 


.20 


-284. 


,05 


610, 


.93 


-315. 


,24 


1.05 


MSk5 [174] 


0. 


,157 


-15. 


,79 


231. 


,17 


385. 


.26 


28. 


,00 


7, 


.57 


-282. 


,55 


607 


.93 


-315. 


,36 


1.05 


MSk5* [119] 


0. 


,156 


-15. 


,78 


243. 


,74 


346. 


,15 


28. 


,00 


7, 


.02 


-290. 


,66 


595, 


.12 


-322. 


,81 


0.80 


MSk6 [174] 


0. 


,157 


-15. 


,79 


231. 


.17 


385. 


.26 


28. 


.00 


9, 


.63 


-274. 


,33 


591, 


.49 


-316. 


,05 


1.05 


MSk7 [175] 


0. 


,157 


-15. 


,80 


231. 


,22 


385. 


,36 


27. 


,95 


9, 


.40 


-274. 


,63 


592 


.08 


-315. 


,38 


1.05 


MSk8 [175] 


0. 


,157 


-15. 


,80 


229. 


.31 


391. 


.01 


27. 


,93 


8, 


.26 


-280. 


,01 


597, 


.59 


-315. 


,49 


1.10 


MSk9 [175] 


0. 


,157 


-15. 


,80 


233. 


,33 


379. 


,16 


28. 


,00 


10. 


,36 


-270. 


,23 


589, 


.06 


-315. 


,57 


1.00 


MSkA [176] 


0. 


,153 


-15. 


,99 


313. 


.33 


138. 


.15 


30. 


,35 


57. 


17 


-135. 


,34 


197, 


.74 


-453. 


13 


0.79 


MSLO [101] 


0. 


,160 


-16. 


,00 


230. 


,00 


380. 


,32 


30. 


,00 


60. 


,00 


-99. 


.33 


224 


.29 


-360. 


11 


0.80 


NRAPR [177] 


0. 


,161 


-15. 


,85 


225. 


,65 


362. 


,54 


32. 


,78 


59. 


,63 


-123. 


,32 


311, 


.61 


-385. 


,32 


0.69 


PRC45 [178] 


0, 


145 


-15. 


,82 


367, 


.58 


-165, 


.69 


51. 


.01 


141. 


,52 


-23. 


.01 


92, 


.05 


-935. 


,89 


1.00 


RATP [179] 


0. 


,160 


-16. 


,05 


239. 


,52 


349. 


,83 


29. 


,26 


32. 


,39 


-191. 


,23 


440, 


.70 


-338. 


,28 


0.67 


Rs [m] 


0, 


158 


-15. 


,59 


237, 


.42 


348, 


.46 


30. 


.82 


86. 


,39 


-9, 


.21 


22, 


.41 


-400. 


,74 


0.78 


SefmOeS [180] 


0. 


,160 


-15. 


,92 


240. 


,11 


347. 


,11 


88. 


,57 


254. 


,43 


-32. 


.10 


59 


.40 


-1190. 


,85 


0.68 


Sefm074 [180] 


0. 


,160 


-15. 


,81 


240. 


.10 


350. 


.15 


33. 


,40 


88. 


,73 


-33. 


.14 


58, 


.41 


-436. 


12 


0.74 


SefmOSl [180] 


0. 


,161 


-15. 


,69 


237. 


.04 


356. 


.66 


30. 


,76 


79. 


,39 


-39. 


.54 


66, 


.74 


-396. 


,41 


0.81 


Sefm09 [180] 


0. 


,161 


-15. 


,55 


240. 


,06 


349. 


,75 


27. 


,78 


69. 


,96 


-40. 


.80 


70, 


.63 


-358. 


,63 


0.90 


Sefml [180] 


0. 


,161 


-15. 


,40 


240. 


.07 


346. 


.34 


24. 


,81 


59. 


,55 


-46. 


.89 


81, 


.53 


-318. 


,28 


1.00 


SGI [181] 


0. 


,154 


-15. 


,89 


261. 


,75 


297. 


,93 


28. 


,33 


63. 


,86 


-51. 


.99 


194, 


.46 


-362. 


,49 


0.61 


SGII [181] 


0. 


,158 


-15. 


,60 


214. 


.65 


380. 


.91 


26. 


,83 


37. 


,63 


-145. 


,90 


330, 


.41 


-304. 


,90 


0.79 


SGOI [182] 


0. 


,168 


-16. 


,63 


361. 


,59 


-37. 


.36 


45. 


,20 


99. 


,76 


-155. 


,64 


144, 


.36 


-764. 


,53 


0.61 


SGOII [182] 


0. 


,168 


-16. 


,70 


253. 


.28 


346. 


.18 


93. 


,98 


246. 


,02 


-119. 


,57 


272, 


.39 


-1259. 


,44 


0.61 


SI [27] 


0. 


155 


-15. 


.99 


370, 


.38 


-152, 


.32 


29. 


,24 


1, 


.22 


-461. 


,84 


141, 


.44 


-469. 


,66 


0.91 


SII [27] 


0. 


,148 


-15. 


,99 


341. 


,40 


-15. 


,76 


34. 


16 


50. 


,02 


-265. 


,72 


104, 


.75 


-568. 


17 


0.58 


SIII [183] 


0. 


,145 


-15. 


,85 


355. 


.37 


-101. 


.38 


28. 


16 


9. 


.91 


-393. 


,73 


130, 


.45 


-456. 


,01 


0.76 


SIII* [184] 


0. 


,148 


-16. 


,07 


361. 


,15 


-107. 


,94 


31. 


,97 


28. 


,70 


-358. 


,37 


84, 


.84 


-539. 


13 


0.78 


SIV [183] 


0. 


,151 


-15. 


.96 


324. 


.55 


68. 


.84 


31. 


,22 


63. 


,50 


-136. 


,72 


79, 


.45 


-504. 


,22 


0.47 


Ski' [185] 


0. 


,155 


-15. 


,99 


370. 


,38 


-152. 


,32 


29. 


,35 


35. 


,34 


-259. 


16 


141, 


.44 


-485. 


,71 


0.91 


SK255 [68] 


0, 


157 


-16. 


,33 


254, 


.93 


350, 


.09 


37. 


,40 


95. 


,05 


-58, 


.33 


94, 


.23 


-498. 


11 


0.80 


SK272 [68] 


0, 


155 


-16. 


,28 


271, 


.51 


305, 


.31 


37. 


,40 


91. 


,67 


-67, 


.78 


134, 


.36 


-514. 


,70 


0.77 


SkA [186] 


0. 


155 


-15. 


,99 


263, 


.16 


300, 


.13 


32. 


,91 


74. 


,62 


-78, 


.46 


174 


.54 


-441. 


,08 


0.61 


Ska25s20 [187] 


0. 


,161 


-16. 


,07 


220. 


.75 


413. 


.45 


33. 


,78 


63. 


,81 


-118. 


,22 


314, 


.13 


-381. 


,56 


0.98 


Ska35sl5 [187] 


0. 


,158 


-16. 


,01 


238. 


,89 


378. 


,88 


30. 


,56 


30. 


,60 


-222. 


,90 


481, 


.99 


-357. 


,96 


1.01 


Ska35s20 [187] 


0. 


,158 


-16. 


,08 


240. 


.27 


378. 


.65 


33. 


,57 


64. 


,83 


-120. 


,32 


284, 


.54 


-407. 


11 


1.00 


Ska35s25 [187] 


0. 


,158 


-16. 


14 


241. 


,30 


378. 


,94 


36. 


,98 


98. 


,89 


-23. 


.57 


97 


.46 


-461. 


,60 


0.99 


Ska45s20 [187] 


0. 


,156 


-16. 


,08 


260. 


,21 


330. 


,55 


33. 


,39 


66. 


,21 


-119. 


,99 


251, 


.77 


-433. 


13 


1.02 


SkB [186] 


0. 


,155 


-15. 


,99 


263. 


.16 


300. 


.13 


23. 


.88 


47. 


,54 


-78. 


.46 


174, 


.54 


-309. 


,50 


0.61 


Skll [188] 


0. 


,160 


-15. 


,95 


242. 


,75 


346. 


,14 


37. 


,53 


161. 


,05 


234. 


,67 


-328 


.02 


-502. 


,01 


0.69 


SkI2 [188] 


0. 


,158 


-15. 


,78 


240. 


.93 


339. 


.70 


33. 


,37 


104. 


,33 


70. 


,69 


51, 


.62 


-408. 


,21 


0.68 


SkI3 [188] 


0. 


,158 


-15. 


,98 


258. 


,19 


303. 


,86 


34. 


,83 


100. 


,53 


73. 


,04 


211 


.54 


-411. 


,80 


0.58 
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TABLE II - continued from previous page 



Model 




Po 


Eo 


K 


O 


K- 


1 


J 




L 




-^sym 


Qsym 






m* 


SkI4 [188| 


0, 


.160 


-15.95 


247, 


.95 


331, 


.21 


29, 


.50 


60. 


,39 


-40, 


.56 


351 


.16 


-322. 


,23 


0.65 


SkI5 [188] 


0, 


.156 


-15.85 


255, 


.79 


301, 


.95 


36, 


.64 


129. 


,33 


159, 


,57 


11 


.73 


-463. 


,74 


0.58 


SkI6 [189| 


0, 


.159 


-15.89 


248, 


.17 


326, 


.58 


29 


.90 


59. 


,24 


-46, 


.77 


378 


.12 


-324. 


,26 


0.64 


SkM [122] 


0, 


.160 


-15.77 


216, 


.61 


386, 


.09 


30, 


.75 


49. 


,34 


-148, 


.81 


323 


.34 


-356. 


,91 


0.79 


SkM* [190] 


0, 


.160 


-15.77 


216, 


.61 


386, 


.09 


30 


.03 


45. 


,78 


-155, 


,94 


330, 


.47 


-349. 


,00 


0.79 


SkMl [191] 


0, 


.160 


-15.77 


216, 


.61 


386, 


.09 


25 


.17 


-35, 


.37 


-388, 


,89 


912, 


.87 


-239. 


,72 


0.79 


SkMP [192] 


0, 


.157 


-15.56 


230, 


.87 


338, 


.05 


29, 


.89 


70, 


,31 


-49, 


,82 


159, 


.44 


-368. 


,73 


0.65 


SkO [193| 


0, 


.160 


-15.84 


223, 


.34 


392, 


.86 


31 


.97 


79, 


14 


-43, 


.17 


131, 


.13 


-378. 


,80 


0.90 


SkO' [193] 


0, 


.160 


-15.75 


222, 


.36 


390, 


.83 


31, 


.95 


68, 


,94 


-78, 


.82 


223, 


.37 


-371. 


,29 


0.90 


SkP [194] 


0, 


.163 


-15.95 


200, 


.97 


435, 


.43 


30 


.00 


19, 


,68 


-266, 


.60 


508, 


.35 


-342. 


,04 


1.00 


SKRA [195] 


0, 


.159 


-15.78 


216, 


.98 


378, 


.76 


31, 


.32 


53, 


,04 


-139, 


,28 


310, 


.84 


-364. 


,92 


0.75 


SkSl [196| 


0, 


.161 


-15.86 


228, 


.43 


382, 


.76 


28 


.75 


30, 


,52 


-218, 


,69 


379, 


.24 


-350. 


,66 


0.86 


SkS2 [196| 


0, 


.161 


-15.89 


229, 


.02 


382, 


.73 


29 


.23 


37, 


,84 


-218, 


,07 


270, 


.03 


-381. 


,86 


0.85 


SkS3 [196| 


0, 


.161 


-15.88 


228, 


.83 


382, 


.62 


28, 


.84 


51, 


.74 


-157, 


.38 


154, 


.06 


-381. 


,30 


0.85 


SkS4 [196] 


0, 


.163 


-15.88 


228, 


.08 


385, 


.45 


28 


.35 


23, 


.28 


-238, 


.42 


438, 


.06 


-338. 


,77 


0.87 


SkSCl [197] 


0, 


.161 


-15.85 


234, 


.58 


380, 


.50 


28, 


.10 


0, 


.13 


-312, 


.03 


673, 


.32 


-312. 


,62 


1.00 


SkSC2 [197] 


0, 


.161 


-15.90 


235, 


.13 


381, 


.60 


24, 


.74 


11, 


,00 


-228, 


,22 


505, 


.69 


-276. 


,35 


1.00 


SkSCS [197] 


0, 


.161 


-15.85 


234, 


.49 


380, 


.32 


27, 


.01 


0, 


.81 


-296, 


,20 


641, 


.65 


-299. 


,75 


1.00 


SkSC4 \W§] 


0, 


.161 


-15.87 


234, 


.72 


380, 


.79 


28, 


.80 


-2, 


.12 


-329, 


,49 


708, 


.23 


-320. 


,20 


1.00 


SkSC4o [199] 


0, 


.161 


-15.87 


234, 


.74 


380, 


.79 


27, 


.00 


-9, 


.67 


-338, 


,03 


725, 


.33 


-295. 


,70 


1.00 


SkSCS [200] 


0, 


.161 


-15.85 


234, 


.50 


380, 


.34 


30, 


.99 


-6, 


.97 


-375, 


,08 


799, 


.41 


-344. 


,58 


1.00 


SkSC6 [200] 


0, 


.161 


-15.92 


235, 


.41 


382, 


.13 


24, 


.57 


11, 


,00 


-226, 


,26 


501, 


.80 


-274. 


,39 


1.00 


SkSClO [200| 


0, 


.161 


-15.96 


235, 


.89 


383, 


.08 


22, 


.83 


19, 


13 


-172, 


,77 


394, 


.81 


-256. 


,47 


1.00 


SkSCll [201] 


0, 


.161 


-15.87 


234, 


.72 


380, 


.79 


28, 


.80 


-2, 


.12 


-329, 


,49 


708, 


.23 


-320. 


,20 


1.00 


SkSC14 [199] 


0, 


.161 


-15.92 


235, 


.41 


382, 


.13 


30, 


.00 


33, 


13 


-202, 


,83 


454, 


.93 


-347. 


,84 


1.00 


SkSC15 [199] 


0, 


.161 


-15.88 


234, 


.93 


381, 


.17 


28, 


.00 


6, 


.72 


-284, 


,47 


618, 


.21 


-313. 


,89 


1.00 


SkSP.l [119] 


0, 


.162 


-15.90 


230, 


.02 


502, 


.64 


28, 


.00 


7, 


.17 


-289, 


,55 


662, 


.66 


-316. 


,92 


0.80 


SkT [202| 


0, 


.148 


-15.40 


333, 


.36 


-29, 


.01 


33, 


.66 


80, 


,83 


-78, 


.93 


69, 


.87 


-570. 


,95 


0.60 


SkTl [113] 


0, 


.161 


-15.98 


236, 


.16 


383, 


.52 


32, 


.02 


56, 


18 


-134, 


,83 


318, 


.99 


-380. 


,68 


1.00 


SkT2 [113] 


0, 


.161 


-15.94 


235, 


.73 


382, 


.67 


32, 


.00 


56, 


16 


-134, 


,67 


318, 


.66 


-380. 


,48 


1.00 


SkT3 [113] 


0, 


.161 


-15.95 


235, 


.74 


382, 


.70 


31, 


.50 


55, 


,31 


-132, 


,05 


313, 


.43 


-374. 


14 


1.00 


SkT4 [113] 


0, 


.159 


-15.96 


235, 


.50 


382, 


.94 


35, 


.24 


93, 


,49 


-24, 


.46 


97, 


.84 


-433. 


,36 


1.00 


SkT5 [113] 


0, 


.164 


-16.00 


201, 


.69 


436, 


.81 


37, 


.00 


98, 


,53 


-24, 


.97 


99, 


.88 


-402. 


,76 


1.00 


SkT6 [113] 


0, 


.161 


-15.96 


235, 


.95 


383, 


.15 


29, 


.97 


30, 


,85 


-211, 


,53 


472, 


.36 


-346. 


,54 


1.00 


SkT7 [113] 


0, 


.161 


-15.94 


235, 


.64 


372, 


.22 


29, 


.52 


31, 


12 


-209, 


,85 


439, 


.35 


-347. 


,42 


0.83 


SkT8 [113] 


0, 


.161 


-15.94 


235, 


.70 


372, 


.37 


29, 


.92 


33, 


,72 


-187, 


,52 


476, 


.25 


-336. 


,59 


0.83 


SkT9 [113] 


0, 


.160 


-15.88 


234, 


.91 


370, 


.97 


29, 


.76 


33, 


,74 


-185, 


,62 


471, 


.98 


-334. 


,76 


0.83 


SkTl* [113] 


0, 


.162 


-16.20 


238, 


.95 


388, 


.75 


32, 


.31 


56, 


,58 


-136, 


,66 


322, 


.86 


-384. 


,07 


1.00 


SkT3* [113] 


0, 


.162 


-16.20 


238, 


.95 


388, 


.76 


31, 


.97 


56, 


,32 


-133, 


,65 


316, 


.82 


-379. 


,93 


1.00 


SkTla [180] 


0, 


.161 


-15.98 


236, 


.16 


383, 


.52 


32, 


.02 


56, 


18 


-134, 


,83 


318, 


.99 


-380. 


,68 


1.00 


SkT2a [180] 


0, 


.161 


-15.94 


235, 


.73 


382, 


.67 


32, 


.00 


56, 


16 


-134, 


,67 


318, 


.66 


-380. 


,48 


1.00 


SkT3a [180] 


0, 


.161 


-15.95 


235, 


.74 


382, 


.70 


31, 


.50 


55, 


,31 


-132, 


,05 


313, 


.43 


-374. 


14 


1.00 


SkT4a [180] 


0, 


.159 


-15.96 


235, 


.50 


382, 


.94 


35, 


.45 


94, 


13 


-24, 


.46 


97, 


.84 


-436. 


,19 


1.00 
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TABLE II - continued from previous page 



Model 




Po 


E 


O 


K 


O 


K' 


1 


J 




L 




-^sym 


Qsyi 


:n 




r 


m* 


SkT5a fl80] 


0. 


,164 


-16, 


,00 


201. 


.69 


436. 


.81 


37. 


,00 


98. 


,53 


-24. 


.97 


99. 


,88 


-402. 


,76 


1.00 


SkT6a [180] 


0. 


,161 


-15, 


,96 


235. 


.95 


383. 


.15 


29. 


,97 


30. 


,85 


-211. 


,53 


472. 


,36 


-346. 


,54 


1.00 


SkT7a [180] 


0. 


,161 


-15. 


,94 


235. 


,64 


372. 


,22 


29. 


,52 


31. 


12 


-209. 


,85 


439. 


,35 


-347. 


,42 


0.83 


SkTSa [180] 


0. 


,161 


-15, 


,94 


235. 


.70 


372. 


.37 


29. 


,92 


33. 


,72 


-187. 


,52 


476. 


,25 


-336. 


,59 


0.83 


SkT9a [180] 


0. 


,160 


-15. 


,88 


234. 


,91 


370. 


,97 


29. 


,76 


33. 


,74 


-185. 


,62 


471. 


,98 


-334. 


,76 


0.83 


SkTK [203] 


0. 


,168 


-16. 


,70 


253. 


,28 


346. 


,18 


35. 


,57 


41. 


,59 


-221. 


,79 


527. 


,94 


-414. 


,46 


0.61 


SKX [204] 


0. 


,155 


-16. 


,05 


271. 


.06 


297. 


.42 


31. 


,10 


33. 


18 


-252. 


12 


379. 


,69 


-414. 


,81 


0.99 


SKXce [204] 


0. 


155 


-15. 


,86 


268, 


.19 


294, 


.59 


30. 


,15 


33. 


,48 


-238. 


,39 


356. 


,93 


-402. 


,51 


1.01 


SKXm [204] 


0. 


,159 


-16. 


.04 


238. 


.09 


380. 


.38 


31. 


,20 


32. 


,08 


-242. 


,76 


428. 


,73 


-384. 


,00 


0.97 


Skxsl5 [205] 


0. 


,161 


-15. 


,76 


201. 


,10 


424. 


,57 


31. 


,88 


34. 


,79 


-197. 


10 


516. 


,30 


-332. 


,38 


0.97 


Skxs20 [205] 


0. 


,162 


-15. 


,81 


201. 


.95 


425. 


.56 


35. 


,50 


67. 


,06 


-122. 


,31 


328. 


,52 


-383. 


,37 


0.96 


Skxs25 [205] 


0. 


,161 


-15. 


,87 


202. 


,92 


-440. 


,88 


39. 


,60 


100. 


10 


-50. 


.28 


145. 


,99 


-440. 


,88 


0.96 


Skz-1 [128] 


0. 


,160 


-16. 


,01 


230. 


.08 


365. 


,25 


32. 


,00 


54. 


14 


-184. 


,08 


217. 


,03 


-422. 


,99 


0.70 


SkzO [128] 


0. 


,160 


-16. 


,01 


230. 


.08 


365. 


.24 


32. 


,00 


35. 


10 


-242. 


,20 


405. 


16 


-397. 


,08 


0.70 


Skzl [128] 


0. 


,160 


-16. 


,01 


230. 


,08 


365. 


,25 


32. 


,01 


27. 


,67 


-242. 


,40 


535. 


,38 


-364. 


,50 


0.70 


Skz2 [128] 


0. 


,160 


-16. 


.01 


230. 


.07 


365. 


.23 


32. 


,01 


16. 


,81 


-259. 


,66 


682. 


,63 


-333. 


,83 


0.70 


Skz3 [128j 


0. 


,160 


-16. 


,01 


230. 


,09 


365. 


,26 


32. 


,01 


12. 


,96 


-241. 


,91 


794. 


,95 


-299. 


,08 


0.70 


Skz4 [128] 


0. 


,160 


-16. 


,01 


230. 


.08 


365. 


.26 


32. 


,01 


5, 


.75 


-240. 


.86 


923. 


.89 


-266. 


,24 


0.70 


SLyO [206] 


0, 


160 


-15. 


,97 


229, 


.66 


364, 


.01 


31. 


,98 


47. 


,11 


-116. 


,23 


508. 


,68 


-324. 


,23 


0.70 


SLyl [206] 


0, 


160 


-15. 


,99 


229, 


.81 


364, 


.35 


31. 


,99 


47. 


,07 


-116. 


,49 


509. 


,36 


-324. 


,27 


0.70 


SLy2 [206] 


0, 


161 


-15. 


,99 


229, 


.92 


364, 


.21 


32. 


,00 


47. 


,46 


-115. 


,13 


506. 


,52 


-324. 


,69 


0.70 


SLy230a [45] 


0. 


,160 


-15. 


,99 


229. 


,89 


364. 


,18 


31. 


,99 


44. 


,32 


-98. 


.22 


602. 


,87 


-293. 


,91 


0.70 


SLy230b [45] 


0. 


,160 


-15. 


,97 


229. 


.91 


363. 


.10 


32. 


,01 


45. 


,97 


-119. 


,72 


521. 


,50 


-322. 


,92 


0.69 


SLy3 [206] 


0. 


160 


-15. 


,94 


229, 


.51 


362, 


.56 


31. 


,97 


45. 


,36 


-121. 


,90 


524. 


,75 


-322. 


,39 


0.70 


SLy4 [20r| 


0, 


160 


-15. 


,97 


229, 


.91 


363 


.11 


32. 


,00 


45. 


,94 


-119. 


,73 


521. 


,53 


-322. 


,83 


0.69 


SLy5 [202 


0, 


161 


-15. 


,99 


229, 


.92 


364, 


.16 


32. 


,01 


48. 


,15 


-112. 


,76 


500. 


,67 


-325. 


,38 


0.70 


SLy6 [202 


0. 


159 


-15. 


,92 


229, 


.86 


360, 


.24 


31. 


,96 


47. 


,45 


-112. 


,71 


510. 


,63 


-323. 


,03 


0.69 


SLy7 [202 


0, 


158 


-15. 


,90 


229, 


.75 


359, 


.22 


31. 


,99 


46. 


,94 


-114. 


,34 


517. 


14 


-322. 


,60 


0.69 


SLy8 [206] 


0. 


160 


-15. 


,97 


229, 


.89 


363, 


.27 


32. 


,00 


47. 


,18 


-115. 


,59 


509. 


,88 


-324. 


,09 


0.70 


SLy9 [206] 


0, 


151 


-15. 


,80 


229, 


.84 


350, 


.42 


31. 


,98 


54. 


,86 


-81, 


.42 


462. 


,35 


-326. 


,92 


0.67 


SLylO [202 


0. 


156 


-15. 


,90 


229, 


.68 


358, 


.32 


31. 


,90 


38. 


,51 


-142. 


,18 


591. 


,23 


-313. 


17 


0.68 


SQMCl [156| 


0. 


,137 


-14. 


,00 


328. 


,76 


-143. 


.78 


29. 


,68 


-6. 


.70 


-504. 


,25 


218. 


,08 


-461. 


10 


0.93 


SQMC2 [156] 


0. 


,140 


-14. 


.29 


330. 


.10 


-121. 


,75 


28. 


,70 


8. 


.67 


-408. 


,41 


145. 


,55 


-463. 


,63 


0.83 


SQMC3 [156] 


0. 


,161 


-15. 


,98 


366. 


.97 


-130. 


.22 


45. 


,78 


91. 


,80 


-210. 


,95 


163. 


,48 


-794. 


,33 


0.82 


SQMC600 [152 


0. 


,174 


-15. 


,74 


217. 


.00 


388. 


.62 


34. 


,38 


46. 


,38 


-215. 


16 


396. 


,85 


-410. 


,40 


0.81 


SQMC650 [152 


0. 


,172 


-15. 


,57 


218. 


,11 


376. 


,75 


33. 


,65 


52. 


,92 


-173. 


15 


349. 


,74 


-399. 


,28 


0.78 


SQMC700 [157] 


0. 


,171 


-15. 


,49 


222. 


,20 


369. 


.94 


33. 


,47 


59. 


,06 


-140. 


,84 


313. 


,84 


-396. 


,85 


0.76 


SQMC750 [152 


0. 


,171 


-15. 


,60 


222. 


.86 


365. 


,83 


33. 


,75 


64. 


,67 


-117. 


,51 


288. 


,41 


-399. 


,38 


0.74 


SSk [52 


0. 


,161 


-16. 


16 


229. 


.31 


375. 


.38 


33. 


,50 


52. 


,78 


-119. 


15 


482. 


,24 


-349. 


,42 


0.72 


SV [183] 


0. 


,155 


-16. 


,05 


305. 


.70 


175. 


.78 


32. 


,82 


96. 


,09 


24. 


17 


48. 


,00 


-497. 


11 


0.38 


SV-bas [115] 


0. 


,160 


-15. 


,91 


233. 


,45 


379. 


,28 


30. 


,00 


32. 


,37 


-221. 


,75 


410. 


,93 


-363. 


,36 


0.90 


SV-min [115] 


0. 


,161 


-15. 


,91 


221. 


.76 


403. 


.08 


30. 


.66 


44. 


,81 


-156. 


,57 


389. 


,56 


-343. 


,99 


0.95 


SVI [183] 


0. 


,143 


-15. 


,76 


363. 


,64 


-153. 


,50 


26. 


,88 


-7. 


.34 


-471. 


,30 


146. 


,04 


-424. 


,18 


0.95 
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TABLE II - continued from previous page 



Model 




Po 


E 


O 


K 


O 


K- 


1 


J 




L 




-^sym 


Qsym 




m* 


SVIl 




0, 


.143 


-15, 


.79 


366, 


.44 


-164, 


.51 


26, 


.96 


-10, 


.16 


-488, 


.90 


149, 


.74 


-423. 


,36 


1.00 


SV-K218 [115] 


0, 


.161 


-15, 


.90 


218, 


.23 


403, 


.15 


30, 


.00 


34. 


,62 


-206, 


.87 


401, 


.58 


-350. 


,65 


0.90 


SV-K226 [115] 


0, 


.160 


-15 


.90 


225, 


.82 


392, 


.14 


30 


.00 


34. 


,09 


-211, 


.92 


401, 


.84 


-357. 


,27 


0.90 


SV-K241 [115] 


0, 


.159 


-15, 


.91 


241, 


.07 


364, 


.54 


30 


.00 


30. 


,95 


-230, 


.77 


416, 


.01 


-369. 


,66 


0.90 


SV-kapOO [115] 


0, 


.160 


-15 


.90 


233, 


.44 


379, 


.15 


30 


.00 


39, 


.44 


-161 


.78 


446, 


.94 


-334. 


,34 


0.90 


SV-kap02 [115] 


0, 


.160 


-15 


.90 


233, 


.44 


379, 


.21 


30 


.00 


35, 


.54 


-193 


.19 


431, 


.91 


-348. 


,69 


0.90 


SV-kap06 [115] 


0, 


.160 


-15, 


.91 


233, 


.45 


379, 


.33 


30, 


.00 


29. 


.33 


-249, 


.75 


388, 


.84 


-378. 


10 


0.90 


SV-mas07 fll5] 


0, 


.160 


-15, 


.89 


233, 


.54 


356, 


.93 


30 


.00 


52. 


15 


-98 


.77 


365 


.68 


-331, 


.96 


0.70 


SV-mas08 fll5] 


0, 


.160 


-15, 


.90 


233, 


.13 


371, 


.28 


30 


.00 


40. 


15 


-172, 


.38 


397, 


.44 


-349. 


,35 


0.80 


SV-maslO [115] 


0, 


.159 


-15, 


.91 


234, 


.33 


383, 


.22 


30 


.00 


28. 


,04 


-252, 


.50 


408, 


.07 


-374. 


,87 


1.00 


SV-s 


iym28 [115] 


0, 


.163 


-16, 


.47 


240, 


.86 


392, 


.55 


28, 


.47 


6, 


.29 


-305, 


.94 


584, 


.47 


-333. 


,41 


0.90 


SV-s 


iym32 [115] 


0, 


.159 


-15, 


.94 


233, 


.81 


380, 


.11 


32 


.00 


57. 


,07 


-148, 


.79 


257, 


.70 


-398. 


,44 


0.90 


SV-s 


iym34 [115] 


0, 


.159 


-15, 


.97 


234, 


.07 


380, 


.82 


34 


.00 


80. 


,95 


-79 


.08 


111, 


.28 


-433. 


,08 


0.90 


SV-tls [115] 


0, 


.160 


-15, 


.89 


233, 


.30 


379, 


.03 


30 


.00 


33. 


,22 


-218, 


.42 


403, 


.90 


-363. 


,79 


0.90 


T [171] 


0, 


.161 


-15 


.93 


235, 


.66 


382, 


.44 


28 


.35 


27. 


.18 


-206 


.76 


462 


.91 


-325. 


,76 


1.00 


Til 


[152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.75 


32, 


.00 


49. 


,46 


-108, 


.76 


486, 


.98 


-326. 


,88 


0.70 


T12 


[152] 


0, 


.161 


-16, 


.00 


230, 


.01 


365, 


.11 


32 


.00 


49. 


,38 


-108, 


.75 


488 


.50 


-326. 


,63 


0.70 


T13 


[152] 


0, 


.161 


-16, 


.00 


230, 


.01 


364, 


.78 


32, 


.00 


49. 


,53 


-108, 


.06 


487 


.57 


-326. 


,69 


0.70 


T14 


[152] 


0, 


.161 


-15, 


.99 


230, 


.01 


364, 


.48 


32, 


.00 


49. 


,48 


-108, 


.12 


488 


.35 


-326. 


,57 


0.70 


T15 


[152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.32 


32 


.00 


49. 


,65 


-107, 


.91 


485 


.83 


-326. 


,95 


0.70 


T16 


[152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.68 


32, 


.00 


49. 


,45 


-108, 


.75 


487 


.24 


-326. 


,83 


0.70 


T21 


[152] 


0, 


.161 


-16, 


.03 


230, 


.01 


366, 


.49 


32 


.00 


49. 


,77 


-108, 


.03 


483 


.25 


-327. 


,37 


0.70 


T22 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


365, 


.95 


32, 


.00 


49. 


,57 


-108, 


.50 


485, 


.74 


-327. 


,04 


0.70 


T23 


[152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.63 


32 


.00 


49. 


,59 


-108, 


.27 


485 


.95 


-326. 


,97 


0.70 


T24 


[152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.37 


32, 


.00 


49. 


,85 


-107, 


.22 


484, 


.00 


-327. 


14 


0.70 


T25 


[152] 


0, 


.161 


-15, 


.99 


230, 


.01 


364, 


.24 


32, 


.00 


49. 


12 


-109, 


.21 


491, 


.85 


-326. 


16 


0.70 


T26 


[152] 


0, 


.161 


-15, 


.98 


230, 


.01 


363, 


.48 


32 


.00 


48. 


,76 


-110, 


.15 


495 


.92 


-325. 


,64 


0.70 


T31 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.28 


32, 


.00 


49. 


,75 


-108, 


.00 


483 


.82 


-327. 


,27 


0.70 


T32 


[152] 


0, 


.161 


-16, 


.03 


230, 


.01 


366, 


.39 


32 


.00 


50. 


,28 


-106, 


.20 


478 


.97 


-327. 


,80 


0.70 


T33 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.10 


32, 


.00 


49. 


,66 


-108, 


.23 


484, 


.88 


-327. 


13 


0.70 


T34 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.28 


32 


.00 


50. 


10 


-106, 


.81 


480 


.71 


-327. 


,60 


0.70 


T35 


[152] 


0, 


.161 


-16, 


.00 


230, 


.01 


364. 


.84 


32, 


.00 


49. 


,59 


-107, 


.85 


487 


.05 


-326. 


,74 


0.70 


T36 


[152] 


0, 


.161 


-15, 


.99 


230, 


.01 


364, 


.51 


32, 


.00 


49. 


,05 


-109, 


.62 


491, 


.98 


-326. 


,20 


0.70 


T41 


[152] 


0, 


.162 


-16, 


.06 


230, 


.01 


368, 


.36 


32 


.00 


50. 


,60 


-106, 


.02 


473 


.67 


-328. 


,60 


0.70 


T42 


[152] 


0, 


.162 


-16, 


.05 


230, 


.01 


368, 


.04 


32, 


.00 


50. 


,70 


-105, 


.51 


473, 


.28 


-328. 


,59 


0.70 


T43 


[152] 


0, 


.162 


-16, 


.04 


230, 


.01 


367, 


.39 


32 


.00 


50. 


,57 


-105, 


.66 


475 


.23 


-328. 


,31 


0.70 


T44 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


365, 


.91 


32, 


.00 


50. 


,05 


-106, 


.76 


481 


.62 


-327. 


,45 


0.70 


T45 


^ 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.10 


32 


.00 


49. 


,66 


-108, 


.24 


484 


.73 


-327. 


16 


0.70 


T46 


[152] 


0, 


.161 


-16, 


.00 


230, 


.01 


364, 


.75 


32 


.00 


49. 


,93 


-106, 


.59 


484 


.25 


-327. 


,00 


0.70 


T51 


[152] 


0, 


.162 


-16, 


.05 


230, 


.01 


367, 


.96 


32, 


.00 


50. 


,69 


-105, 


.52 


473 


.55 


-328. 


,55 


0.70 


T52 


[152] 


0, 


.161 


-16, 


.06 


230, 


.01 


368, 


.07 


32 


.00 


50. 


,68 


-105, 


.55 


473 


.55 


-328. 


,55 


0.70 


T53 


[152] 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.21 


32, 


.00 


50. 


,03 


-106, 


.99 


481, 


.50 


-327. 


,50 


0.70 


T54 


[152] 


0, 


.161 


-16, 


.03 


230, 


.01 


366, 


.73 


32 


.00 


50. 


,27 


-106, 


.36 


478 


.71 


-327. 


,85 


0.70 
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TABLE II - continued from previous page 



Model 




Po 


E 


O 


K 


O 


K' 


1 


J 




L 




-^sym 


Qsym 




m* 


T55 [152] 


0, 


.161 


-16, 


.03 


230, 


.01 


366, 


.66 


32. 


.00 


50. 


,24 


-106. 


,49 


479, 


.02 


-327, 


.83 


0.70 


T56 [152] 


0, 


.161 


-16, 


.01 


230, 


.01 


365, 


.26 


32. 


.00 


50. 


13 


-106. 


19 


481 


.83 


-327, 


.34 


0.70 


T61 [152] 


0, 


.162 


-16, 


.07 


230, 


.01 


368, 


.76 


32, 


.00 


50. 


,79 


-105. 


,56 


471 


.67 


-328, 


.85 


0.71 


T62 [152] 


0, 


.162 


-16, 


.07 


230, 


.01 


368, 


.93 


32. 


.00 


50. 


,33 


-107. 


,25 


475, 


.46 


-328, 


.49 


0.71 


T63 [152] 


0, 


.162 


-16, 


.06 


230, 


.01 


368, 


.30 


32, 


.00 


51. 


,07 


-104. 


,36 


469 


.72 


-329, 


.00 


0.70 


T64 [152] 


0, 


.162 


-16, 


.03 


230, 


.01 


366, 


.74 


32. 


.00 


50. 


,49 


-105. 


,65 


476 


.73 


-328, 


.08 


0.70 


T65 [152] 


0, 


.162 


-16, 


.04 


230, 


.01 


367, 


.37 


32. 


.00 


50. 


,50 


-105. 


,90 


475, 


.82 


-328, 


.25 


0.70 


T66 [152] 


0, 


.161 


-16, 


.02 


230, 


.01 


366, 


.04 


32, 


.00 


50. 


,30 


-105. 


,96 


479 


.25 


-327, 


.72 


0.70 


v070 [208| 


0, 


.157 


-15, 


.78 


230, 


.99 


384, 


.93 


27. 


.98 


-3, 


.45 


-361. 


.15 


591, 


.72 


-346, 


.20 


1.05 


v075 [208| 


0, 


.157 


-15 


.80 


231, 


.29 


385, 


.51 


28, 


.00 


-0, 


.31 


-341, 


.88 


587, 


.67 


-340, 


.52 


1.05 


v080 [208| 


0, 


.157 


-15, 


.79 


231, 


.17 


385, 


.26 


28. 


.00 


2, 


.23 


-325. 


.61 


585, 


.53 


-335, 


.29 


1.05 


v090 [208] 


0, 


.157 


-15 


.79 


231, 


.17 


385, 


.26 


28, 


.00 


5, 


.04 


-304, 


.26 


593 


.46 


-326, 


.10 


1.05 


vlOO [208] 


0, 


.157 


-15 


.79 


231, 


.17 


385, 


.26 


28, 


.00 


8, 


.73 


-281, 


.39 


588 


.25 


-319, 


.22 


1.05 


vl05 [208| 


0, 


.157 


-15, 


.79 


231, 


.17 


385, 


.26 


28. 


.00 


7, 


.08 


-284. 


.51 


611 


.85 


-315, 


.20 


1.05 


vllO [208] 


0, 


.157 


-15 


.79 


231, 


.17 


385, 


.26 


28, 


.00 


7, 


.51 


-279, 


.62 


617 


.86 


-312, 


.19 


1.05 


z [mj 


0, 


.159 


-15, 


.97 


330, 


.30 


64, 


.98 


26. 


,82 


-49, 


.70 


-657. 


.85 


495, 


.24 


-369, 


.43 


0.84 


ZRla [178] 


0, 


.173 


-16 


.99 


398, 


.74 


-186, 


.01 


9. 


,84 


-57, 


.61 


-471. 


,08 


103, 


.34 


-98 


.56 


1.00 


ZRlb [178] 


0, 


.173 


-16, 


.99 


398, 


.74 


-186, 


.01 


18, 


.50 


-31, 


.62 


-471, 


.08 


103, 


.34 


-266, 


.61 


1.00 


ZRlc [178] 


0, 


.173 


-16, 


.99 


398, 


.74 


-186, 


.01 


31. 


,50 


7, 


.36 


-471. 


,08 


103, 


.34 


-518, 


.70 


1.00 


ZR2a [178] 


0, 


.173 


-16 


.99 


324, 


.78 


184, 


.94 


1, 


.62 


-82, 


.36 


-397, 


.29 


474, 


.98 


49, 


.99 


1.00 


ZR2b [rre] 


0, 


.173 


-16, 


.99 


324, 


.78 


184, 


.94 


11. 


.95 


-51, 


.39 


-397. 


.29 


474, 


.98 


-118, 


.22 


1.00 


ZR2c [178] 


0, 


.173 


-16, 


.99 


324, 


.78 


184, 


.94 


27, 


.43 


-4, 


.93 


-397, 


.29 


474, 


.98 


-370, 


.53 


1.00 


ZR3a [178] 


0, 


.175 


-16, 


.99 


198, 


.79 


475, 


.65 


-138. 


.96 


-504, 


.42 


-271, 


.89 


768, 


.04 


1547, 


.69 


1.00 


ZR3b [178] 


0, 


.175 


-16 


.99 


198, 


.79 


475, 


.65 


-100, 


.46 


-388, 


.91 


-271, 


.89 


768, 


.04 


1131, 


.04 


1.00 


ZRSc [178] 


0, 


.175 


-16, 


.99 


198, 


.79 


475, 


.65 


-42. 


,71 


-215, 


.66 


-271. 


.89 


768, 


.04 


506, 


.06 


1.00 


Zs [171] 


0, 


.163 


-15, 


.88 


233, 


.33 


368, 


.95 


26. 


.69 


-29, 


.38 


-401. 


.43 


883, 


.05 


-271, 


.61 


0.78 


Zs* [171] 


0, 


.162 


-15 


.96 


234, 


.87 


369, 


.16 


28, 


.80 


-4, 


.53 


-332, 


.64 


725 


.10 


-312, 


.58 


0.77 
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TABLE III: Constraints for the Skyrme Models. The t symbol on the minus sign indicates that the model satisfies all the constraints except one. 



Model 


SMI 


SM2 SM3 


SM4 PNMl 


PNM2 


MIXl 


MIX2 


MIX3 


MIX4 MIX5 


Model 


SMI 


SM2 


SM3 


SM4 PNMl 


PNM2 


MIXl 


MIX2 


MIX3 


MIX4 


MIX5 


BSkl 


+ 


+ + 


+ 


+ 










+ 


- 


GSl 


+ 


+ 




+ + 


+ 




+ 




+ 


+ 


BSk2 


+ 


+ + 


+ 


+ 










+ 


- 


GS2 




+ 




+ 










+ 




BSk2' 


+ 


+ + 


+ 


+ 










+ 


- 


GS3 








+ 








+ 


+ 




BSk3 


+ 


+ + 


+ 


+ 










+ 


- 


GS4 


+ 


+ 




+ 












+ 


BSk4 


+ 


+ + 


+ 


+ 










+ 


- 


GS5 




+ 




+ 












+ 


BSk5 


+ 


+ + 


+ 


+ 










+ 


- 


GS6 








+ 














BSk6 


+ 


+ + 


+ 


+ 










+ 


+ 


GSkI 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


-t 


+ 


+ 


BSk7 


+ 


+ + 


+ 


+ 










+ 


+ 


GSkll 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


_t 


+ 


+ 


BSkS 


+ 


+ + 


+ 


+ 










+ 


+ 


KDE 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


-t 


+ 


+ 


BSk9 


+ 


+ + 


+ 


+ 




+ 






+ 


+ 


KDEOv 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


.t 


+ 


+ 


BSklO 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


KDEOvl 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


.t 


+ 


+ 


BSkll 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


LNS 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


+ 


+ 


+ 


BSkl2 


+ 


+ 


+ 






+ 






+ 


+ 


MSkl 


+ 


+ 


+ 


+ + 




+ 






+ 


+ 


BSkl3 


+ 


+ 


+ 






+ 






+ 


+ 


MSk2 


+ 


+ 


+ 


+ + 




+ 






+ 


+ 


BSkl4 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


MSk3 


+ 


+ 


+ 


+ + 










+ 




BSkl5 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


MSk4 


+ 


+ 


+ 


+ + 










+ 




BSkl6 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


MSk5 


+ 


+ 


+ 


+ + 










+ 




BSkl7 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


MSk5* 


+ 


+ 




+ 










+ 




BSkl8 


+ 


+ 


+ 


+ 




+ 






+ 


+ 


MSk6 


+ 


+ 


+ 


+ + 










+ 




BSkl9 


+ 


+ 


+ 


+ 




+ 






+ 


- 


MSk7 


+ 


+ 


+ 


+ + 










+ 




BSk20 


+ 


+ 


+ 


+ 


+ 


+ 






+ 


_ 


MSk8 


+ 


+ 


+ 


+ + 










+ 




BSk21 




-4- 


+ 


+ 






4- 




+ 


+ 


MSk9 


-4- 


4- 


4- 


4- 4- 










4- 




E 
















+ 






MSkA 








+ 




+ 


+ 


+ 


+ 


+ 


Es 


+ 


+ 


+ 
















MSLO 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


.t 


+ 


+ 


f- 


+ 


+ + 


+ 




+ 


+ 


+ 




+ 


+ 


NRAPR 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


+ 


+ 


+ 


f+ 


+ 


+ + 


+ 




+ 


+ 


+ 




+ 




PRC45 




















+ 


fO 


+ 


+ + 


+ 




+ 


+ 


+ 




+ 


+ 


RATP 


+ 


+ 




+ + 










+ 




FPLyon 


+ 


+ + 


+ 


+ 




+ 


+ 




+ 


+ 


Rs 


+ 


+ 




+ 


+ 


+ 




+ 




+ 


Gs 


+ 


+ 


+ 




+ 


+ 




+ 




+ 


Sefm068 


+ 


+ 




+ 












+ 
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TABLE III - continued from previous page 



A/TnHpl 


SMI 


SM2 SMS 


SM4 PNMl 


PNM2 


MIXl 


MIX2 


MTXS MTX4 MTX'i 

xvx±,/v*_i ±vxxyv*± ±V±±,/V^J 


MnHpl 


SMI 


SM2 


SM3 


SM4 PNMl 


PNM2 


MIXl 


MIX2 


MIX3 


MTX4 MIX'S 


OcilliU / ^ 




1 


1 




1 




1 1 


SkM 


1 


1 


1 


1 1 




1 


1 




1 1 




_i_ 


_l_ 


_l_ 




_|_ 




4- - -U 


SkM* 


_|_ 


_|_ 


_|_ 


_|_ _|_ 




_l_ 


_|_ 




_|_ _|_ 




_L 


-I- 










- - 4- 


SkMl 








4- 












Sefml 


_L 












_ 4_ 4- 


SkMP 


4- 


4- 


4- 








4- 










1 


1 1 






1 


1 1 




1 


1 


1 


1 




1 




1 


1 


SGII 


_l_ 


_1_ _L 


_|_ _|_ 








_|_ 


SkO' 


_l_ 


_|_ 


_|_ 


_l_ _|_ 


.t 


_l_ 


_|_ 


_i_ 


_l_ _|_ 
















1 1 




1 


1 


1 


1 1 




1 






1 


SGOII 


+ 


+ 


+ 


- 


- 


- 


+ 


SKRA 


+ 


+ 


+ 


+ + 


+ 


+ 


+ 


.t 


+ + 


SI 














4- 4- - 


SkSl 


4- 


4- 


4- 


4- 4- 










_1_ 


SII 














_i_ 4_ _ 
-f- -f- - 


SkS2 


_|_ 


_|_ 


_|_ 










_i_ 


_1_ 


SIII 






-U 








4- 4- - 


SkS3 


_|_ 


_|_ 


_|_ 


_|_ 






_|_ 


_i_ 


_l_ 


SIII* 






_L 








_L. _|_ _1_ 

_l_ _l_ _j_ 


SkS4 


_i_ 


_|_ 


_|_ 


_i_ _i_ 










4- 


SIV 












_i_ 


4- 4- - 


SkSCl 


4- 




4- 


4- 










4- 


QTV-ocrr 
OlvZOO 


1 

-r 


1 

_|_ 


00 

1 1 
-r "T 


1 
1 

"T 






1 1 1 
_l_ _^ 

1 1 




1 
1 

-r 


1 
1 

-r 


1 
1 

-r 


1 1 

1 

-r 










1 
1 


iDlvZ ( Z 




1 


1 1 


1 






1 1 




1 


1 


1 












1 






1 


1 1 




1 


1 


1 1 1 
_l_ 




1 


1 


1 


1 




















4- 






_L. _L 4_ 




4- 






4- 




4- 


















1 




1 1 


OKOv_^0 


1 


1 


1 


1 1 










1 




1 




1 1 


1 


1 




1 1 1 




1 


1 


1 


1 1 










1 1 




1 


1 


1 1 


1 






1 1 




1 


1 


1 


1 










1 






1 


1 1 




1 


1 


1 1 1 
T 




1 


1 


1 


1 1 




1 






1 1 


SkB 






4- 






4- 


- 4- - 








_i_ 












4- 


Skll 


+ 


+ 


+ 








+ - + 


SkSP.l 


+ 


+ 


+ 


+ + 










+ 


SkI2 


+ 


+ 


+ 




+ 




+ - + 


SkT 








+ 


+ 


+ 




+ 


+ 


SkI3 


+ 


+ 


+ + 




+ 




+ - + 


SkTl 


+ 


+ 


_t 


+ + 


+ 


+ 


+ 


+ 


+ + 


SkI4 


+ 


+ 


+ + 






+ 


+ + 


SkT2 


+ 


+ 


.t 


+ + 


+ 


+ 


+ 


+ 


+ + 


SkI5 


+ 


+ 


+ 








+ - + 


SkT3 


+ 


+ 


-t 


+ + 


+ 


+ 


+ 


+ 


+ + 


SkI6 


+ 


+ 


+ + 






+ 


+ + 


SkT4 


+ 


+ 


+ 


+ + 


+ 






+ 


+ 


Continued on next page 



TABLE III - continued from previous page 



TVTnHpl 


SMI SM2 SM3 


SM4 PNMl PNM2 


MIXl MIX2 


MIX3 MIX4 MIX5 


MoHpI 


SMI 


SM2 SM3 SM4 PNMl 


PNM2 


MIXl 


MIX2 MIX3 


MIX4 MIX5 


SkT5 


-1- -1- -k 


4- - 4- 




+ 


+ 


SLyO 


4_ 


4-4-4- 4- 


4- 


4- 


1 _t 


4- 4- 


SkT6 


_i_ _i_ 


4-4- 




+ 


+ 


ST.vl 
oijy ± 


4_ 


4-4-4- 4- 


4_ 


4_ 




4_ 4_ 


SkT7 


_L _U _ 
_^ _ 


4-4- 




+ 




SLy2 


4_ 


4_ 4_ 4_ 4_ 
_j_ _l_ _^ 


4- 


4_ 


4- -t 


4_ 4_ 


SkT8 


_L 4- _ 


4-4- 




+ 


+ 


\Di-iy ^OKJa 


4- 


4_ 4_ 4_ 
_l_ 


4_ 


4_ 


4- 


i_ 4_ 


SkT9 


-1- A- 


4_ _i_ 

1 




+ 


+ 




4_ 


4-4-4- 4- 
1^ 


4_ 


4_ 


1 t 


4_ 4_ 

1 


SkTl* 


+ + -t 


4- -I- 4- 


4_ 4_ 


+ + 


+ 


SLy3 


4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


SkT3* 


+ + -t 


4-4- 4- 
\ \ \ 


4_ 4_ 


+ + 


+ 




4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


SkTla 


+ + -t 


+ + + 


4_ 4_ 


+ + 


+ 


SLy5 


4_ 


4_ 4_ 4^ 4_ 


4- 


4_ 


+ -t 


4_ 4_ 




+ + -t 


4-4- 4- 

\ \ \ 


4_ 4_ 


+ + 


+ 


SLy6 


4- 


4-4-4- 4- 


4_ 


4- 


+ -t 


4_ 4_ 


SkT3a 


+ + -+ 


+ + + 


4_ 4^ 


+ + 


+ 




4- 


-1- -1- -1- -1- 

\ \ \ \ 


4- 


4- 


+ _t 


4- 4- 


SkT4a 


-I- -I- -l- 


4-4- 4- 




+ 


+ 


SLyS 


4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


SkTSa 


-1- -1- -1- 

\ \ \ 


-1- - -1- 




+ 


+ 


SLvQ 


4- 


-1- -1- -1- -1- 

\ \ \ \ 


4- 


4- 


+ -t 


4- 4- 




_L _L _ 

\ \ 


4-4- 




+ 


+ 


SLylO 


4- 


4-4-4- 
_l_ 


4- 


4- 




4_ 4_ 


SkT7a 


-4- 4- - 


4- 4-rn 
CO 

1 1 
+ -r 




+ 
















4- 


OK ± Oa 


1 1 
_|_ _|_ - 




+ 


+ 












1 


1 


OK _L y<i 


1 1 


1 1 




+ 


+ 








1 






1 1 


SkTK 


-1- 4- - 


4- _ 


4_ 


+ + 






4- 


4_ 4_ 4_ 4_ 
_l_ _^ _l_ 




4_ 


4- 4- 


4- 


SKX 




4-4- 


4- 


+ + 






4- 


4_ 4_ 4_ 4_ 
_l_ 


4- 


4- 


4- 4- 


4_ 4_ 


9KXrp 
oxvyvcc 


- 4- - 


4-4- 


4- 


+ + 






4_ 


4-4-4- 4- 


4_ 


4_ 


4_ 4_ 


4_ 4_ 


SKXm 


4_ 4_ 


4-4- 


4- 


+ + 






4_ 


4- - 4- 4- 




4_ 


4- 4- 


4_ 4_ 


OK_X.D ±ij 


4-4-4- 


4-4- 


4_ 


+ 


+ 


SSk 


4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


OKa-D^U 


4-4-4- 


4-4- 4- 
_l_ _^ 




+ + 


+ 


SV 




_ _ _ 1 




4- 






OK-X-O^iJ 


4-4-4- 
_l_ _^ _j_ 


4- - 




+ 


+ 


O V UdiO 


4- 


4-4-4- 4- 
_l_ _l_ _j_ _|_ 




4_ 




4- 


Skz-1 


+ + + 


+ + 


+ + 


+ + 




SV-min 


+ 


+ + + + 




+ 


+ 


+ + 


SkzO 


+ + + 


+ + 


+ 


+ + 




SVI 










+ 


+ 


Skzl 


+ + + 


+ + 


+ 


+ 




SVII 










+ 


+ 


Skz2 


+ + + 


+ 


+ 


+ 




SV-K218 


+ 


+ + + + 




+ 




+ 


Skz3 


+ + + 




+ 


+ 




SV-K226 


+ 


+ + + + 




+ 




+ 


Skz4 


+ + + 


+ 


+ 




+ 


SV-K241 


+ 


+ - + + 




+ 




+ 
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TABLE III - continued from previous page 



A/TnHpl 


SMI 


SM2 


SM3 


SM4 PNMl 


PNM2 


MIXl 


MIX2 


MIX3 MIX4 MIX5 


IVTnHpl 


SMI 


SM2 SlVfi SM4 PNMl 


PNM2 


MIXl 


MTX2 MTXS 


MTX4 MIX'S 


SV-kanOn 


-l_ 


-1- 


-\- 


+ 


+ 




-(_ 




+ 


+ 


T41 


-(_ 


4-4-4- 4- 

\ 


-|- 


-|- 


+ _t 


+ + 


SV-kann2 


_i_ 


_|_ 


_j_ 


+ 


+ 




_i_ 




+ 


+ 


T42 


4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


RV-kanflR 


_L. 


_L. 


_L. 


+ 


+ 




_L. 




+ + 




T43 




_L. _L J_ 4. 
_^ _l_ _l_ 


4- 


_L 


+ -t 


_^ 


O V -IIld.bU / 


1 


1 




+ 


+ 


1 


1 


1 


+ 


+ 


T44 


1 


4-4-4- 4- 


1 


1 


4- t 


1 1 


iD V -illcL&Uo 


1 


1 


1 


+ 


+ 




1 


1 


+ 


+ 


-L 


1 


4- 4- 4_ 4- 


1 


1 




1 1 


iD V -lll(lo_LU 


1 


1 


1 


+ 


+ 




1 




+ + 




T46 


1 


4- 4- 4_ 4- 


1 


1 


-k t 


1 1 


o V -by iiizo 


1 


1 




+ 


+ 








+ 




± o± 


1 


4- 4- 4_ 4- 


1 


1 


-k t 


1 1 


SV-sym32 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ + 


+ 


T52 


+ 


+ + + + 


+ 


+ 


+ -t 


+ + 


o V -byiiio'Tt 


1 


1 




+ 






1 




+ + 


+ 


-L OO 


1 


4-4-4- -U 


1 


1 


1 t 


1 1 


SV-tls 


_l_ 

I 


4- 

V 


_i_ 


+ 


+ 








+ 




T54 




4-4-4- 4- 


_|_ 




1 t 


_l_ _|_ 


T 


_l_ 
1 


_l_ 
1 




+ 


+ 








+ 


+ 


T55 


_l_ 


4-4-4- 4- 










Til 


_|_ 


_|_ 


_i_ 


+ 


+ 




_l_ 




-t + 


+ 


T56 


_i_ 


4-4-4- 4- 




_l_ 




_l_ _|_ 


T12 








+ 


+ 


4- 






-t + 


+ 


T61 


4- 


_L. 4_ _l_ 

_^ _l_ _|_ 


4- 


4- 


4- -t 




T1 

± LO 

T14 


1 

_|_ 


1 


1 


+ 
+ 


+ 0) 

o 

+ 


1 

_L 
r 


1 

_l_ 


1 

t 

1 


-t + 

-t + 


+ 
+ 


T62 

T63 


1 

_l_ 


4- 4- -J- 4- 

_^ _l_ _^ 

4-4-4- 4- 


1 

_l_ 


1 

_l_ 


4- t 
-\- t 


1 1 
4- 4- 


T15 


_|_ 


_l_ 


_l_ 


+ 


+ 


_|_ 


_l_ 


_l_ 


-t + 


+ 


T64 


4_ 


4-4-4- 4- 


4_ 


4_ 


+ -t 


4_ 4_ 


T16 




4_ 




+ 


+ 


_i_ 


_1_ 


_}_ 


-t + 


+ 


T65 


4- 


_L. J_ _L. _l_ 

_l_ _l_ _^ _l_ 


4- 


-L 


+ -t 




T21 


1 


1 


1 


+ 


+ 


1 


1 


1 


-t + 


+ 


J. uo 


1 


4- 4- -J- -J- 

_l_ _|_ 


1 


1 


4- t 
-|- 


1 1 


T22 




1 


1 


+ 


+ 


1 


1 


1 


-t + 


+ 


VU ( u 


1 


4-4-4- 4- 








1 


T23 


_L 

r 




_i_ 


+ 


+ 


_|_ 


_l_ 




-t + 


+ 


v075 




4-4-4- 4- 








_|_ 


T24 




_l_ 


_l_ 


+ 


+ 


_|_ 


_l_ 


_l_ 


-t + 


+ 


v080 




4-4-4- 4- 








_|_ 


T25 






_L 


+ 


+ 






_L 


-t + 


+ 


v090 


4- 


_J_ J_ _L. -J- 

-f- -f- ~f- -f- 








4- 


T26 


1 


1 


1 


+ 


+ 


1 


1 


1 


-t + 


+ 


V -LUU 


1 


4- _U -1- -J- 

_^ _l_ _|_ 








1 


T31 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


-t + 


+ 


vl05 


+ 


+ + + + 








+ 


T32 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


-t + 


+ 


vllO 


+ 


+ + + + 








+ 


T33 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


-t + 


+ 


Z 












+ 


T34 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


-t + 


+ 


ZRla 




+ 








+ 


T35 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


-t + 


+ 


ZRlb 




+ 








+ 


T36 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


-t + 


+ 


ZRlc 








+ 


+ 


+ 
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TABLE III - continued from previous page 



Model SMI SM2 SMS SM4 PNMl PNM2 MIXl MIX2 MIX3 MIX4 MIX5 Model SMI SM2 SMS SM4 PNMl PNM2 MIXl MIX2 MIXS MIX4 MIX5 
ZR2a _ _ _ _ + - _ _ _ _ + ZRSb - + + + - - - - - - + 

ZR2b _ - - - + - _ _ - - + ZRSc - + + + - - - - - - + 

ZR2c ----- - - - - + - Zs + + - + - - ----- 

ZRSa - + + + - - - - - - + Zs* + + - + - 



TABLE IV: List of the Skyrme parameterizations and their compliance with macroscopic constraints used in this work. For SMI, SM2 and MIXl-5 

we give the standard deviation; the constraint is satisfied if the standard deviation is less of equal one. For SMS, SM4, PNMl and PNM2 we give the 
fraction of full density range in % in which the constraint is not satisfied. Letter L (H) indicates that the Skyrme prediction starts to fail in the beginning 
(end) of the density range. No indication of letters means that the model fails in the middle of the density range. For models that fail in 100%, the 
letter U (D) indicates that its curve is above (below) the band defined by the constraint. The t symbol marks a constraint which is the only one which 
is not satisfied (out of 11) by a particular Skyrme parameterization. For more details see text. 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 


MIX3 


MIX4 MIX5 


BSkl 


0.04 


-0.63 











lOOD 


-1.88 


-2.82 


1.30 


0.69 


1.70 


BSk2 


0.12 


-0.64 











lOOD 


-1.80 


-2.78 


1.21 


0.69 


3.80 


BSk2' 


0.11 


-0.64 











lOOD 


-1.80 


-2.79 


1.21 


0.69 


3.90 


BSkS 


0.16 


-0.64 











lOOD 


-1.83 


-2.85 


1.18 


0.69 


4.40 


BSk4 


0.23 


-0.67 


2.0L 








lOOD 


-1.80 


-2.53 


1.26 


0.54 


2.00 


BSk5 


0.24 


-0.66 


2.7L 








lOOD 


-1.52 


-2.03 


1.19 


0.31 


1.20 


BSk6 


-0.03 


-0.66 











lOOD 


-1.80 


-2.29 


1.43 


0.54 


0.90 


BSk7 


-0.02 


-0.66 











lOOD 


-1.80 


-2.22 


1.43 


0.54 


0.80 


BSkS 


0.01 


-0.66 








2.4 


lOOD 


-1.80 


-2.40 


1.44 


0.62 


0.90 


BSk9 


0.04 


-0.65 











67.8L 


-1.00 


-1.09 


1.26 


0.00 


0.30 


BSklO 


0.29 


-0.66 


9.6L 








lOOD 


-1.00 


-1.15 


1.06 


-0.08 


0.80 


BSkll 


0.27 


-0.66 


7.0L 








lOOD 


-1.00 


-1.09 


1.06 


-0.08 


0.70 


BSkl2 


0.27 


-0.66 


6.9L 








lOOD 


-1.00 


-1.11 


1.06 


-0.08 


0.80 


BSkl3 


0.27 


-0.66 


7.0L 








lOOD 


-1.00 


-1.07 


1.06 


-0.15 


0.70 


BSkl4 


0.31 


-0.68 


11.4L 








0.0 


-1.00 


-0.78 


1.12 


-0.23 


0.60 


BSkl5 


0.39 


-0.67 


14.6L 








lOOD 


-1.00 


-1.36 


1.14 


0.08 


0.80 


BSkl6 


0.39 


-0.67 


14.0L 








lOOD 


-1.00 


-1.28 


1.14 


0.00 


0.70 


BSkl7 


0.39 


-0.67 


14.0L 








lOOD 


-1.00 


-1.21 


1.14 


0.00 


0.70 


BSkl8 


0.39 


-0.67 


14.1L 








lOOD 


-1.00 


-1.21 


1.15 


0.00 


0.70 


BSkig 


0.24 


-0.80 


23.6L 








lOOD 


-1.00 


-1.45 


1.15 


0.15 


1.90 


BSk20 


0.38 


-0.84 


30.6L 








0.0 


-1.00 


-1.15 


1.28 


0.08 


1.30 


BSk21 


0.53 


-0.85 


33. 4L 








17.1L 


-1.00 


-0.64 


1.55 


0.00 


0.70 


E 


3.45 


-1.27 


lOOU 


lOOU 


lOOU 


lOOD 


-1.94 


-4.96 


0.91 


1.46 


-1.20 


Es 


0.62 


-0.70 


29.2L 





lOOU 


lOOD 


-2.42 


-5.27 


1.43 


2.00 


-1.20 


f_ 


n no 
u.uu 


-0.59 





n 
u 


1 AT /II 9TT 


n 
u 


-u.zu 


-u. 1 y 


1 42 


0.00 


0.10 


f+ 


0.00 


-0.59 








19.0L/25.8H 





-0.20 


-0.91 


1.40 


0.08 


1.10 


fO 


0.00 


-0.59 








16.8L/18.8H 





-0.20 


-0.87 


1.41 


0.08 


0.70 


FPLyon 


-0.43 


-0.60 











8.2L 


-0.63 


-0.85 


1.30 


-0.08 


0.20 


Gs 


0.24 


-0.70 


7.8L 





90.2H 





-0.55 


1.96 


0.81 


-1.46 


0.40 


GSl 


0.17 


0.22 


92.3H 


0.7H 








-1.46 


-0.43 


1.96 


-0.38 


0.20 


GS2 


2.34 


-0.76 


50.5L 


43.6L/16.6 





11.8L 


-2.62 


-1.54 


1.18 


-0.08 


1.30 


GS3 


5.66 


-2.26 


lOOU 


lOOU 





lOOD 


-4.40 


-3.25 


0.93 


0.46 


93.00 


GS4 


0.17 


0.29 


lOOD 


10.9H 





lOOD 


-7.87 


-4.26 


2.32 


2.23 


1.00 


GS5 


2.31 


-0.68 


40.9L 


31.4L/5.2 





lOOD 


-5.52 


-3.90 


1.72 


1.23 


-0.90 


GS6 


5.70 


-2.17 


lOOU 


lOOU 





lOOD 


-7.27 


-5.61 


1.92 


2.54 


-1.30 


GSkI 


0.01 


-0.59 











0.7L 


-0.19 


0.30 


1.04t 


-0.62 


0.20 


GSkll 


0.11 


-0.60 











4.0 


-0.80 


-0.52 


1.03+ 


-0.38 


0.90 



Continued on next page 



62 



TABLE IV - continued from previous page 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 


MIX3 


MIX4 MIX5 


KDE 


-0.20 


-0.64 














-0.21 


-0.92 


1.27+ 


0.00 


0.30 




-0.04 


-0.65 





Q 


Q 


Q 


0.19 


-0.71 


1.15+ 


-0.08 


0.40 


IX 1 I I J\f V X 


-0.08 


-0.63 














0.83 


-0.18 


1.05+ 


-0.23 


0.30 


LNS 


-0.64 


-0.63 














0.37 


0.19 


0.94 


-0.54 


0.60 


MSkl 


0.12 


-0.64 


0.0 








lOOD 


-1.00 


-1.34 


1.12 


0.00 


0.40 


MSk2 


0.06 


-0.63 


0.0 








lOOD 


-1.00 


-1.37 


1.12 


0.08 


0.40 




0.11 


-0.64 


0.0 








lOOD 


-1.80 


-2.83 


1.30 


0.77 


2.00 


MSk4 


0.04 


-0.63 


0.0 








lOOD 


-1.80 


-2.82 


1.29 


0.69 


1.70 


MSk5 


0.04 


-0.63 


0.0 








lOOD 


-1.80 


-2.80 


1.29 


0.69 


1.60 




0.46 


-0.71 


14.9L 





6.4H 


lOOD 


-1.80 


-2.83 


1.25 


0.69 


5.40 


MSkfi 


0.04 


-0.63 


0.0 








lOOD 


-1.80 


-2.69 


1.29 


0.69 


1.30 


MSk7 


0.04 


-0.63 


0.0 








lOOD 


-1.82 


-2.70 


1.29 


0.69 


1.30 


MSk8 


-0.02 


-0.62 











lOOD 


-1.83 


-2.76 


1.29 


0.69 


1.30 


MSk9 


0.11 


-0.64 











lOOD 


-1.80 


-2.65 


1.29 


0.62 


1.30 


MSkA 


2.78 


-1.12 


lOOU 


lOOU 





48. 4L 


-0.86 


-0.05 


0.58 


-0.62 


0.70 


MST.n 


0.00 


-0.64 











Q 


-1.00 


0.11 


1.06+ 


-0.69 


0.60 


NR APT! 


-0.14 


-0.67 














0.11 


0.09 


0.93 


-0.54 


0.80 


PRC45 


4.59 


-1.73 


lOOU 


lOOU 


26. OH 


lOOU 


7.40 


4.64 


-1.91 


-1.31 


0.10 


RATP 


0.32 


-0.70 


16.3L 








lOOD 


-1.30 


-1.42 


1.17 


0.08 


1.70 


Rs 


0.25 


-0.70 


8.4L 





86. 4H 





-0.67 


1.58 


0.85 


-1.31 


0.50 




0.34 


-0.71 


17.8L 


0.0 


88. 6H 


lOOU 


22.43 


10.91 


-3.22 


-1.38 


0.20 


Sefmn74 


0.34 


-0.70 


17.2L 


0.0 


65.4 


0.0 


0.36 


1.71 


0.67 


-1.23 


0.50 




0.23 


-0.69 


12.6L 


0.0 


80.2 


15. 4L 


-0.70 


1.19 


0.87 


-1.15 


0.50 


Scfm09 


0.34 


-0.70 


17.3L 


0.0 


80. 2H 


59. 4L 


-1.89 


0.66 


1.07 


-1.08 


0.50 


Scfml 


0.34 


-0.71 


17.8L 


0.0 


80. 8H 


38. 8L 


-3.08 


0.09 


1.28 


-1.00 


0.50 


SGI 


1.06 


-0.80 


37.6L 








52. 6L 


-1.67 


0.33 


1.05 


-0.85 


0.70 


SGII 


-0.51 


-0.64 











lOOD 


-2.27 


-1.13 


1.35 


-0.15 


1.10 


SGOI 


4.39 


-1.47 


lOOU 


lOOU 


10.0L/46.4H 


lOOU 


5.08 


2.32 


-1.02 


-0.92 


0.70 


SGOII 


0.78 


-0.71 


41. 7L 


0.0 


lOOU 


lOOU 


24.59 


10.45 


-3.57 


-1.23 


0.30 


SI 


4.68 


-1.70 


lOOU 


lOOU 


10 0/71 2H 

l.\J.\J/ § ±.^±1. 


lOOD 


-1.30 


-3.15 


0.50 


0.54 


25.20 


SII 


3.71 


-1.43 


lOOU 


lOOU 








0.66 


-0.44 


-0.01 


-0.46 


1.20 


SIII 


4.18 


-1.60 


lOOU 


lOOU 





lOOD 


-1.74 


-2.67 


0.57 


0.31 


4.10 


SIII* 


4.37 


-1.62 


lOOU 


lOOU 


0.0 


lOOD 


-0.21 


-1.63 


0.14 


-0.08 


1.00 


SIV 


3.15 


-1.26 


lOOU 


lOOU 





54. 7L 


-0.51 


0.31 


0.32 


-0.77 


1.30 


Ski' 


4.68 


-1.70 


lOOU 


lOOU 








-1.26 


-1.26 


0.41 


-0.23 


0.90 




0.83 


-u. / u 










0.0 


1.96 


2.06 


0.35 


1 1 


u.ou 


SK272 


1.38 


-0.79 


45. IL 








0.0 


1.96 


1.87 


0.26 


-1.08 


0.50 


SkA 


1.11 


-0.80 


39.6L 








42.7L 


0.16 


0.92 


0.64 


-0.92 


0.80 


Ska25s20 


-0.31 


-0.57 














0.51 


0.32 


0.95 


-0.54 


0.20 


Ska35sl5 


0.30 


-0.64 


5.1L 








lOOD 


-.78 


-1.52 


1.07 


0.15 


0.40 


Ska35s20 


0.34 


-0.64 


7.6Lt 











0.43 


0.38 


0.82 


-0.62 


0.20 


Ska35s25 


0.38 


-0.64 


8.5L 











1.79 


2.27 


0.54 


-1.23 


0.10 


Ska45s20 


1.01 


-0.74 


33.6L 








12.5L 


0.36 


0.46 


0.68 


-0.69 


0.20 
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TABLE IV - continued from previous page 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 MIX3 


MIX4 MIX5 


SkB 


1.11 


-0.80 


39.6L 





26.8 


24. IL 


-3.45 


-0.58 


1.32 


-0.69 


1.30 


Skll 


0.42 


-0.71 


21.2L 


Q 


98. 4H 


lOOU 


2.01 


5.73 


0.33 


-2.46 


0.10 


SkI2 


0.36 


-0.72 


15.3L 





80.4H 


71 OT./I 5 5H 


0.35 


2.57 


0.81 


-1.46 


0.10 


SkI3 


0.94 


-0.79 


38.3L 








lOOU 


0.93 


2.36 


0.79 


-1.23 


0.00 


SkI4 


0.60 


-0.74 


28.5L 








43. 9L 


-1.20 


0.13 


1.26 


-0.54 


0.10 


SkI5 


0.86 


-0.80 


33.9L 





71.8 


lOOU 


1.66 


3.96 


0.53 


-1.77 


0.00 


SkI6 


0.61 


-0.75 


27.4L 





0.0 


55. 3L 


-1.04 


0.07 


1.25 


-0.54 


0.10 


SkM 


-0.45 


-0.63 











5.6 


-0.70 


-0.48 


1.08 


-0.38 


0.90 


SkM* 


-0.45 


-0.63 











lOOD 


-0.99 


-0.68 


1.12 


-0.31 


1.00 


SkMl 


-0.45 


-0.63 


0.0 





lOOU 


lOOD 


-2.93 


-5.19 


1.68 


2.15 


-1.30 


SkMP 


0.03 


-0.72 








40.0 


55. 2L 


-1.04 


0.68 


1.02 


-0.92 


0.70 


SkO 


-0.22 


-0.61 








64.0 


57. IL 


-0.21 


1.17 


0.96 


-1.08 


0.30 


SkO' 


-0.25 


-0.62 











23.6Lt 


-0.22 


0.61 


1.00 


-0.77 


0.30 


SkP 


-0.97 


-0.53 











lOOD 


-1.00 


-2.13 


1.15 


0.38 


2.30 


SKRA 


-0.43 


-0.64 














-0.47 


-0.28 


1.04t 


-0.38 


0.90 


SkSl 


-0.05 


-0.63 





Q 


Q 


lOOD 


-1.50 


-1.53 


1.11 


0.00 


1.80 


SkS2 


-0.03 


-0.63 











lOOD 


-1.31 


-1.12 


0.95 


-0.23 


2.00 


SkS3 


-0.04 


-0.63 








54.2 


56.4H 


-1.46 


-0.35 


0.95 


-0.62 


1.40 


SkS4 


-0.06 


-0.63 








0.0 


lOOD 


-1.66 


-1.93 


1.17 


0.23 


2.20 


SkSCl 


0.15 


-0.64 


2.5L 





20 nT,/67 SH 


lOOD 


-1.76 


-3.21 


1.31 


0.92 


5.00 




0.17 


-0.64 


3.4L 


Q 


0.0 


lOOD 


-3.10 


-2.61 


1.49 


0.62 


1.30 




0.15 


-0.64 


2.3L 





10. 8H 


lOOD 


-2.20 


-3.18 


1.37 


0.92 


17.40 


SkSC4 


0.16 


-0.64 


2.7L 





lOOU 


lOOD 


-1.48 


-3.34 


1.27 


1.00 


-6.70 




0.16 


-0.64 


3.0L 





lOOU 


lOOD 


-2.20 


-3.76 


1.39 


1.23 


-1.50 




0.15 


-0.64 


2.3L 





lOOU 


lOOD 


-0.60 


-3.61 


1.14 


1.08 


-2.00 




0.18 


-0.64 


4.1L 


Q 


Q 


lOOD 


-3.17 


-2.61 


1.50 


0.62 


1.30 




0.20 


-0.63 


4.9L 








lOOD 


-3.87 


-2.16 


1.60 


0.31 


0.70 


SkSCl 1 


0.16 


-0.64 


2.7L 





lOOU 


lOOD 


-1.48 


-3.34 


1.27 


1.00 


-6.70 


0-IA.0 V_ _L rt 


0.18 


-0.64 


4.1L 








lOOD 


-1.00 


-1.38 


1.12 


0.08 


0.40 


SkSCl R 


0.16 


-0.64 


3.3L 





2. OH 


lOOD 


-1.80 


-2.85 


1.30 


0.77 


2.10 


SkSP.l 


0.00 


-0.39 











lOOD 


-1.80 


-2.82 


1.28 


0.69 


12.20 


SkT 


3.45 


-1.46 


lOOU 


lOOU 








0.46 


1.27 


-0.03 


-1.08 


0.60 


SkTl 


0.21 


-0.63 


e.OLt 











-0.19 


-0.10 


0.96 


-0.46 


0.30 


SkT2 


0.19 


-0.63 


5.4Lt 











-0.20 


-0.10 


0.96 


-0.46 


0.30 


SkT3 


0.19 


-0.63 


5.3Lt 











-0.40 


-0.15 


0.99 


-0.46 


0.30 


SkT4 


0.18 


-u.oo 


n 
u 











1.10 


1 Q7 


u.oo 


1 1 


n 1 n 


SkT5 


-0.94 


-0.53 








90.2 





1.80 


2.25 


0.84 


-1.23 


0.10 


SkT6 


0.20 


-0.63 


5.3L 








lOOD 


-1.01 


-1.51 


1.13 


0.08 


0.50 


SkT7 


0.19 


-0.66 


7.0L 








lOOD 


-1.19 


-1.49 


1.13 


0.08 


1.20 


SkT8 


0.19 


-0.66 


7.2L 








lOOD 


-1.03 


-1.35 


1.18 


0.08 


0.60 


SkT9 


0.16 


-0.66 


5.1L 








lOOD 


-1.10 


-1.35 


1.19 


0.08 


0.60 


SkTl* 


0.30 


-0.62 


11.7Lt 











-0.08 


-0.08 


0.94 


-0.46 


0.30 


SkT3* 


0.30 


-0.62 


11.7Lt 











-0.21 


-0.09 


0.96 


-0.46 


0.30 
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TABLE IV - continued from previous page 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 


MIX3 


MIX4 MIX5 


SkTla 


0.21 


-0.63 


6.0Lt 











-0.19 


-0.10 


0.96 


-0.46 


0.30 




0.19 


-0.63 


5.4Lt 


Q 


Q 


Q 


-0.20 


-0.10 


0.96 


-0.46 


0.30 


kJJV -L (JCl 


0.19 


-0.63 


5.3Lt 











-.40 


-0.15 


0.99 


-0.46 


0.30 


SkT4a 


0.18 


-0.63 














1.18 


2.01 


0.67 


-1.23 


0.10 


SkT'ia 


-0.94 


-0.53 








90.2 





1.80 


2.25 


0.84 


-1.23 


0.10 




0.20 


-0.63 


5.3L 





Q 


lOOD 


-1.01 


-1.51 


1.13 


0.08 


0.50 


SkT7a 


0.19 


-0.66 


7.0L 








lOOD 


-1.19 


-1.49 


1.13 


0.08 


1.20 


SkT8a 


0.19 


-0.66 


7.2L 








lOOD 


-1.03 


-1.35 


1.18 


0.08 


0.60 


SkTQa 


0.16 


-0.66 


5.1L 








lOOD 


-1.10 


-1.35 


1.19 


0.08 


0.60 


SkTK 


0.78 


-0.71 


41. 7L 


0.0 


10 0L/3S 8H 


76. 7L 


1.23 


-0.91 


0.78 


0.00 


1.60 


SKX 


1.37 


-0.81 


45.4L 








lOOD 


-0.56 


-1.38 


0.78 


0.00 


1.30 


SKXre 


1.27 


-0.81 


43.3L 








lOOD 


-0.94 


-1.36 


0.84 


-0.08 


1.30 


SKXm 


0.27 


-0.64 


5.3L 








lOOD 


-0.52 


-1.44 


0.94 


0.08 


1.30 


SkxslR 


-0.96 


-0.55 











lOOD 


-.25 


-2.29 


1.20 


0.15 


0.50 




-0.94 


-0.55 














1.20^ 


0.50 


0.94 


-0.54 


0.30 




-0.90 


-0.55 








30.2 


5 31. /3 7 


2.84 


2.34 


0.64 


-1.08 


0.20 


Skz-1 


0.00 


-0.67 











65. IH 


-0.20 


-0.21 


0.74 


-0.54 


1.60 


SkzO 


0.00 


-0.67 











lOOD 


-0.20 


-1.27 


0.87 


0.00 


2.30 


Skzl 


0.00 


-0.67 











lOOD 


-0.20 


-1.68 


1.04 


0.23 


2.10 


Skz2 


0.00 


-0.67 








lOOU 


lOOD 


-0.20 


-2.29 


1.20 


0.62 


2.50 


Skz3 


0.00 


-0.67 








lOOU 


lOOD 


-0.20 


-2.50 


1.38 


0.77 


1.30 


Skz4 


0.00 


-0.67 








lOOU 


lOOD 


-0.20 


-2.90 


1.55 


1.08 


-0.70 


SLyO 


-0.01 


-0.67 














-0.21 


-0.60 


1.25t 


-0.08 


0.30 


ST,v1 


-0.01 


-0.67 














-0.20 


-0.61 


1.25t 


-0.08 


0.30 


SLy2 


0.00 


-0.67 














-0.20 


-0.59 


1.24t 


-0.08 


0.30 




0.00 


-0.67 








lOOU 





-0.20 


-0.76 


1.40 


0.08 


-0.20 




0.00 


-0.67 














-0.20 


-0.67 


1.25t 


-0.08 


0.30 


SLy3 


-0.02 


-0.67 














-0.21 


-0.70 


1.26t 


-0.08 


0.30 




0.00 


-0.67 














-0.20 


-0.67 


1.25+ 


-0.08 


0.30 


SLy5 


0.00 


-0.67 














-0.20 


-0.55 


1.24t 


-0.15 


0.30 


SLy6 


0.00 


-0.68 














-0.22 


-0.59 


1.25t 


-0.08 


0.30 


SLv7 


-0.01 


-0.68 














-0.20 


-0.61 


1.25t 


-0.08 


0.30 


SLy8 


0.00 


-0.67 














-0.20 


-0.60 


1.25t 


-0.08 


0.30 




-0.01 


-0.70 














-0.21 


-0.17 


1.23+ 


-0.31 


0.20 


SLylO 


-0.01 


-0.68 








18.8L/29.6H 


4.6L 


-0.24 


-1.08 


1.30 


0.15 


0.30 


SOMCl 


3.29 






lOOU 


95. gH 


lOOD 


-1.13 


-3.59 


0.54 


u.uy 


7 90 


SQMC2 


3.34 


-1.64 


lOOU 


lOOU 


30.8H 


lOOD 


-1.52 


-2.74 


0.53 


0.38 


4.40 


SQMC3 


4.57 


-1.66 


lOOU 


lOOU 


lOOU 


0.0 


5.31 


1.88 


-1.18 


-0.77 


0.50 


SQMC600 


-0.43 


-0.62 


0.0 








lOOD 


0.75 


-0.65 


0.80 


-0.15 


1.30 


SQMC650 


-0.40 


-0.65 


0.0 











0.46 


-0.28 


0.86 


-0.31 


1.00 


SQMC700 


-0.26 


-0.66 


4.1L 











0.39 


0.06 


0.87 


-0.46 


0.90 


SQMC750 


-0.24 


-0.67 


7.9L 








16.0L 


0.50 


0.37 


0.86 


-0.62 


0.70 


SSk 


-0.02 


-0.65 














0.40 


-0.29 


1.12+ 


-0.23 


0.30 
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TABLE IV - continued from previous page 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 MIX3 


MIX4 MIX5 


SV 


2.52 


-1.05 


lOOU 


lOOU 





lOOU 


0.13 


2.12 


0.36 


-1.38 


1.10 




0.11 


-0.64 











lOOD 


-1.00 


-1.42 


1.04 


0.00 


1.30 




-0.39 


-0.59 











lOOD 


-1.00 


-1.30 


1.11 


0.00 


1.30 


SV-K22fi 


-0.14 


-0.62 











lOOD 


-1.00 


-1.33 


1.08 


0.00 


1.30 


SV-K241 


0.37 


-0.67 


13. OL 








lOOD 


-1.00 


-1.50 


1.01 


0.08 


1.40 


SV-kanOn 


0.11 


-0.64 











99. 7L 


-1.00 


-1.03 


1.19 


-0.08 


0.20 


SV-kap02 


0.11 


-0.64 











lOOD 


-1.00 


-1.25 


1.12 


0.00 


0.70 


SV-kap06 


0.11 


-0.64 











lOOD 


-1.00 


-1.59 


0.97 


0.08 


2.10 


SV-mas07 


0.12 


-0.69 


5.4L 











-1.00 


-0.33 


1.21 


-0.38 


0.50 




0.10 


-0.66 











lOOD 


-1.00 


-0.99 


1.12 


-0.15 


0.90 


kJ V L1.10jO±\J 


0.14 


-0.63 











lOOD 


-1.00 


-1.66 


0.99 


0.08 


1.70 


SV-min 


-0.27 


-0.59 











lOOD 


-0.74 


-0.73 


1.14 


-0.23 


0.40 


RV-svm28 

kJ V iSj'llliiO 


0.36 


-0.61 


16.2L 








lOOD 


-1.61 


-2.87 


1.20 


0.69 


7.10 


kJ V Oy ILLtJ^ 


0.13 


-0.64 














-0.20 


-0.05 


0.86 


-0.54 


0.80 


9V-qvm'^4 


0.14 


-0.64 








17.6 


45. 5L 


0.60 


1.28 


0.69 


-1.00 


0.50 


SV-tIs 


0.11 


-0.64 











lOOD 


-1.00 


-1.38 


1.04 


0.00 


1.30 


SVI 


4.45 


-1.71 


lOOU 


lOOU 


27.4H 


lOOD 


-2.25 


-3.63 


0.73 


0.77 


-4.00 


SVII 


4.55 


-1.73 


lOOU 


lOOU 


41. 2H 


lOOD 


-2.22 


-3.79 


0.74 


0.85 


-2.70 


T 


0.19 


-0.64 


5.9L 








lOOD 


-1.66 


-1.71 


1.24 


0.15 


0.50 


Til 


0.00 


-0.67 














-0.20 


-0.47 


1.23+ 


-0.15 


0.30 


T12 


0.00 


-0.67 














-0.20 


-0.48 


1.23t 


-0.15 


0.30 


T13 


0.00 


-0.67 














-0.20 


-0.47 


1.23t 


-0.15 


0.30 


T14 


0.00 


-0.67 














-0.20 


-0.47 


1.23t 


-0.15 


0.30 


T15 


0.00 


-0.67 














-0.20 


-0.46 


1.23t 


-0.15 


0.30 


T16 


0.00 


-0.67 














-0.20 


-0.47 


1.23t 


-0.15 


0.30 


T21 


0.00 


-0.67 














-0.20 


-0.46 


1.23t 


-0.15 


0.30 


T22 


0.00 


-0.67 














-0.20 


-0.47 


1.23t 


-0.15 


0.30 


T23 


0.00 


-0.67 














-0.20 


-0.47 


1.23t 


-0.15 


0.30 


T24 


0.00 


-0.67 














-0.20 


-0.45 


1.23+ 


-0.15 


0.30 


T25 


0.00 


-0.67 








Q 





-0.20 


-0.49 


1.24t 


-0.15 


0.30 


T26 


0.00 


-0.67 














-0.20 


-0.51 


1.24t 


-0.15 


0.30 


T31 


0.00 


-0.67 














-0.20 


-0.46 


1.23t 


-0.15 


0.30 


T32 


0.00 


-0.67 














-0.20 


-0.43 


1.23t 


-0.15 


0.30 


T33 


0.00 


-0.67 














-0.20 


-0.46 


1.23+ 


-0.15 


0.30 


T34 


0.00 


-0.67 














-0.20 


-0.44 


1.23+ 


-0.15 


0.30 


T35 


0.00 


-U.O 1 


u 











-0.20 




1 

J..ZO ' 


-U. J-O 


u.ou 


T36 


0.00 


-0.67 














-0.20 


-0.50 


1.24t 


-0.15 


0.30 


T41 


0.00 


-0.66 


0.5L 











-0.20 


-0.41 


1.22+ 


-0.23 


0.30 


T42 


0.00 


-0.66 


O.IL 











-0.20 


-0.41 


1.22+ 


-0.23 


0.30 


T43 


0.00 


-0.67 














-0.20 


-0.41 


1.23+ 


-0.23 


0.30 


T44 


0.00 


-0.67 














-0.20 


-0.44 


1.23+ 


-0.15 


0.30 


T45 


0.00 


-0.67 














-0.20 


-0.46 


1.23+ 


-0.15 


0.30 


T46 


0.00 


-0.67 














-0.20 


-0.45 


1.23+ 


-0.15 


0.30 
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TABLE IV - continued from previous page 



Model 


SMI 


SM2 SM3(%) 


SM4(%) 


PNM1(%) 


PNM2(%) 


MIXl 


MIX2 


MIX3 


MIX4 MIX5 


T51 


0.00 


-0.66 














-0.20 


-0.41 


1.22+ 


-0.23 


0.30 


T52 


0.00 


-0.66 














-0.20 


-0.41 


1.22+ 


-0.23 


0.30 


T53 


0.00 


-0.67 














-0.20 


-0.44 


1.23+ 


-0.15 


0.30 


T54 


0.00 


-0.67 














-0.20 


-0.43 


1.23+ 


-0.15 


0.30 


T55 


0.00 


-0.67 














-0.20 


-0.43 


1.23+ 


-0.15 


0.30 


T56 


0.00 


-0.67 














-0.20 


-0.44 


1.23+ 


-0.15 


0.30 


T61 


0.00 


-0.66 


0.3L 











-0.20 


-0.40 


1.22+ 


-0.23 


0.30 


T62 


0.00 


-0.66 


0.4L 











-0.20 


-0.43 


1.22+ 


-0.15 


0.30 


T63 


0.00 


-0.66 


0.3L 











-0.20 


-0.39 


1.22+ 


-0.23 


0.30 


T64 


0.00 


-0.67 














-0.20 


-0.42 


1.23+ 


-0.15 


0.30 


T65 


0.00 


-0.67 














-0.20 


-0.42 


1.23+ 


-0.15 


0.30 


T66 


0.00 


-0.67 














-0.20 


-0.43 


1.23+ 


-0.15 


0.30 


v070 


0.03 


-0.63 








4.0H 


lOOD 


-1.81 


-3.41 


1.13 


0.92 


15.50 


v075 


0.04 


-0.63 











lOOD 


-1.80 


-3.24 


1.16 


0.85 


44.70 


v080 


0.04 


-0.63 











lOOD 


-1.80 


-3.10 


1.19 


0.77 


17.10 


v090 


0.04 


-0.63 











lOOD 


-1.80 


-2.94 


1.24 


0.77 


5.20 


vlOO 


0.04 


-0.63 











lOOD 


-1.80 


-2.74 


1.27 


0.69 


1.90 


vl05 


0.04 


-0.63 











lOOD 


-1.80 


-2.83 


1.29 


0.69 


1.70 


vllO 


0.04 


-0.63 











lOOD 


-1.80 


-2.81 


1.31 


0.69 


1.10 


z 


3.34 


-1.27 


lOOU 


lOOU 


lOOU 


lOOD 


-2.27 


-5.98 


1.01 


1.92 


-0.70 


ZRla 


5.62 


-1.77 


lOOU 


lOOU 


Q 


lOOD 


-9.06 


-6.42 


2.41 


4.62 


-0.30 


ZRlb 


5.62 


-1.77 


lOOU 


lOOU 


2.4L 


lOOD 


-5.60 


-4.98 


1.54 


1.69 


-0.50 


ZRlc 


5.62 


-1.77 


lOOU 


lOOU 


lOOU 


lOOD 


-.40 


-2.81 


0.24 


0.38 


2.00 


ZR2a 


3.16 


- -L.UO 




lOOU 





lOOD 


12.35 


-7.80 


3.18 


oy.uo 




ZR2b 


3.16 


-1.03 


lOOU 


lOOU 





lOOD 


-8.22 


-6.08 


2.31 


4.00 


-0.30 


ZR2c 


3.16 


-1.03 


lOOU 


lOOU 


lOOU 


lOOD 


-2.03 


-3.50 


1.01 


0.85 


-3.00 


ZR3a 


-1.04 


-0.45 








lOOD 


lOOD 


68.58 


31.25 


10.90 


-2.08 


0.00 


ZR3b 


-1.04 


-0.45 








99.6H 


lOOD 


53.18 


24.83 


8.75 


-2.31 


0.00 


ZR3c 


-1.04 


-0.45 








91. 2H 


lOOD 


30.08 


15.20 


5.53 


-3.31 


-0.10 


Zs 


0.11 


-0.66 


9.0L 





lOOU 


lOOD 


-2.32 


-4.85 


1.52 


1.85 


-1.50 


Zs* 


0.16 


-0.66 


10.5L 





lOOU 


lOOD 


-1.48 


-3.47 


1.31 


1.08 


10.10 
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TABLE V: Number of Skyrme models consistent with individual macroscopic constraints. 



Constraints SMI SM2 SM3 SM4 PNMl PNM2 MIXl MIX2 MIX3 MIX4 MIX5 
No CSkP 202 214 145 211 186 94 139 102 65 198 159 



TABLE VI: Parameters of the Skyrme interactions consistent with the macroscopic constraints, to 
in MeV.fm^, ti, t2 in MeV.fm^; in McV.fm^+'^'^\ xq, xi, X2, x^, and ai are dimensionless. For 



all parameterizations t4 = X4 = and = x^ = 0. 



Skyrme 


*o 


ti 


t2 






t33 


XQ 




X2 


3:31 


3:32 


3:33 


o-l 


0-2 


0-3 


GSkI 


-1855.5 


397.2 


264.6 


13858.0 


-2694.1 


-319.9 


0.12 


-1.76 


-1.81 


0.13 


-1.19 


-0.46 


0.33 


0.67 


1.00 


GSkll 


-1856.0 


393.1 


266.1 


13842.9 


-2689.7 




0.09 


-0.72 


-1.84 


-0.10 


-0.35 




0.33 


0.67 




KDEOvl 


-2553.1 


411.7 


-419.9 


14603.6 






0.65 


-0.35 


-0.93 


0.95 






0.17 






LNS 


-2485.0 


266.7 


-337.1 


14588.2 






0.06 


0.66 


-0.95 


-0.03 






0.17 






MSLO 


-2118.1 


395.2 


-64.0 


12875.7 






-0.07 


-0.33 


1.36 


-0.23 






0.24 






NRAPR 


-2719.7 


417.6 


-66.7 


15042.0 






0.16 


-0.05 


0.03 


0.14 






0.14 






Ska25s20 


-2180.5 


281.5 


-160.4 


14577.8 






0.14 


-0.80 


0.00 


0.06 






0.25 






Ska35s20 


-1768.8 


263.9 


-158.3 


12904.8 






0.13 


-0.80 


0.00 


0.01 






0.35 






SKRA 


-2895.4 


405.5 


-89.1 


16660.0 






0.08 


0.00 


0.20 


0.00 






0.14 






SkTl 


-1794.0 


298.0 


-298.0 


12812.0 






0.15 


-0.50 


-0.50 


0.09 






0.33 






SkT2 


-1791.6 


300.0 


-300.0 


12792.0 






0.15 


-0.50 


-0.50 


0.09 






0.33 






SkT3 


-1791.8 


298.5 


-99.5 


12794.0 






0.14 


-1.00 


1.00 


0.08 






0.33 






Skxs20 


-2885.2 


302.7 


-323.4 


18237.5 






0.14 


-0.26 


-0.61 


0.05 






0.17 






SQMC650 


-2462.7 


436.1 


-151.9 


14154.5 






0.13 


0.00 


0.00 


0.00 






0.17 






SQMC700 


-2429.1 


371.0 


-96.7 


13773.6 






0.10 


0.00 


0.00 


0.00 






0.16 






SV-sym32 


-1883.3 


319.2 


197.3 


12559.5 






0.01 


-0.59 


-2.17 


-0.31 






0.30 
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TABLE VII: Properties of nuclear matter at saturation as calculated using the Skyrme parameteri- 
zations consistent with the macroscopic constraints. All entries are in MeV except for dimensionless 



Model 




Po 


E 


o 


K 


O 


K' 


J 


L 




S(Po/2) 
J 


3PpNM 


m* 


GSkI 


0. 


,159 


-16, 


.02 


230 


.21 


405, 


.58 


32, 


.03 


63, 


.45 


-364. 


,19 


0. 


,62 


1, 


.02 


0, 


.78 


GSkll 


0. 


,159 


-16 


.12 


233 


.40 


398, 


.73 


30 


.49 


48, 


.63 


-366. 


,54 


0. 


,65 


1, 


.09 





.79 


KDEOvl 


0. 


,165 


-16 


.23 


227, 


.54 


384, 


.86 


34 


.58 


54, 


.69 


-362. 


.78 


0. 


,67 


1, 


.03 





.74 


LNS 


0. 


,175 


-15, 


.32 


210, 


.78 


382, 


.55 


33, 


.43 


61, 


,45 


-384, 


,55 


0, 


,63 


1, 


.06 


0, 


.83 


MSLO 


0. 


,160 


-16, 


.00 


230, 


.00 


380, 


.32 


30, 


.00 


60, 


.00 


-360, 


,11 


0, 


,61 


1, 


.06 


0, 


.80 


NRAPR 


0. 


,161 


-15, 


.85 


225, 


.65 


362 


.54 


32 


.78 


59, 


.63 


-385. 


,32 


0. 


,63 


1, 


.08 


0, 


.69 


Ska25s20 


0. 


161 


-16, 


.07 


220, 


.75 


413 


.45 


33 


.78 


63, 


.81 


-381. 


,56 


0. 


,63 


1, 


.02 





.98 


Ska35s20 


0. 


,158 


-16 


.08 


240, 


.27 


378, 


.65 


33 


.57 


64, 


.83 


-407. 


,11 


0. 


,62 


1, 


.02 


1, 


.00 


SKRA 


0. 


159 


-15, 


.78 


216, 


.98 


378, 


.76 


31, 


.32 


53, 


,04 


-364, 


,92 


0, 


,65 


1, 


.09 





.75 


SkTl 


0. 


,161 


-15, 


.98 


236, 


.16 


383, 


.52 


32, 


.02 


56, 


.18 


-380. 


,68 


0. 


,64 


1, 


.03 


1, 


.00 


SkT2 


0. 


161 


-15, 


.94 


235, 


.73 


382, 


.67 


32, 


.00 


56, 


.16 


-380. 


.48 


0. 


.64 


1, 


.03 


1, 


.00 


SkT3 


0. 


161 


-15, 


.95 


235, 


.74 


382, 


.70 


31, 


.50 


55, 


.31 


-374. 


.14 


0. 


,64 


1, 


.03 


1, 


.00 


Skxs20 


0. 


162 


-15, 


.81 


201, 


.95 


425, 


.56 


35, 


.50 


67. 


.06 


-383. 


.37 


0. 


.63 


1, 


.03 


0, 


.96 


SQMC650 


0. 


,172 


-15, 


.57 


218, 


.11 


376, 


.75 


33 


.65 


52, 


,92 


-399, 


,28 


0, 


,66 


1, 


.10 


0, 


.78 


SQMC700 


0. 


,171 


-15, 


.49 


222, 


.20 


369, 


.94 


33 


.47 


59, 


.06 


-396, 


,85 


0. 


,64 


1, 


.09 


0, 


.76 


SV-sym32 


0. 


159 


-15, 


.94 


233, 


.81 


380, 


.11 


32 


.00 


57, 


.07 


-398. 


,44 


0. 


,63 


1, 


.08 


0, 


.90 
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TABLE VIII: Selected properties on nuclear matter as predicted by consistent Skyrme parameteri- 
zations at 3po- L{po) is included for a comparison. All quantities are in MeV except for m*, which 



is dimensionless. 



Force 


m" 


5 


LiPo) 


L(3po) 




Qsym 


GSkI 


0.536 


60.92 


63.45 


89.71 


-253.92 


947.56 


GSkll 


0.557 


36.42 


48.63 


-69.49 


-731.84 


1226.76 


KDEOvl 


0.493 


55.31 


54.69 


53.66 


-192.75 


1366.50 


LNS 


0.612 


54.15 


61.45 


23.06 


-515.88 


1027.81 


MSLO 


0.571 


53.91 


60.00 


50.80 


-424.29 


776.89 


NRAPR 


0.430 


53.38 


59.63 


27.72 


-468.17 


1024.00 


Ska25s20 


0.942 


58.41 


63.81 


51.87 


-411.53 


1047.02 


Ska35s20 


1.000 


57.75 


64.83 


43.51 


-473.55 


1027.22 


SKRA 


0.748 


44.07 


53.04 


-22.54 


-596.94 


1093.69 


SkTl 


1.000 


48.17 


56.18 


-0.31 


-527.96 


1158.54 


SkT2 


1.000 


48.17 


56.16 


-0.14 


-527.25 


1157.11 


SkT3 


1.000 


47.53 


55.31 


0.78 


-515.94 


1134.50 


Skxs20 


0.899 


61.70 


67.06 


56.83 


-428.23 


1058.16 


SQMC650 


0.779 


33.65 


52.92 


-76.28 


-173.15 


349.74 


SQMC700 


0.507 


49.98 


59.06 


-6.52 


-603.34 


1109.31 


SV-sym32 


0.749 


44.14 


57.07 


-42.43 


-759.99 


1060.52 
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TABLE IX: Relevant combinations of Skyrme parameters of consistent Skyrme parameterizations: 
C£, and Cf are given in MeV.fm^; Q', and Cf' in MeV.fm^; 6^, 6;, 6' and 0'^ in MeV.fm^ 



45|. For full explanation, see text. 



Force 


cS 




Co" 


CI 






©' 

^sym 




GSkI 


-279.2 


68.2 


37.6 


19.3 


602.3 


147.2 


-763.6 


455.28 


GSkll 


-276.9 


154.3 


34.7 


-33.6 


554.6 


546.4 


530.1 


8.17 


KDEOvl 


-282.5 


81.6 


43.3 


14.5 


692.3 


229.8 


-695.0 


462.4 


LNS 


-250.3 


137.9 


25.0 


-19.5 


400.8 


356.4 


267.6 


44.4 


MSLO 


-272.7 


132.7 


32.4 


-23.1 


518.4 


444.3 


296.2 


74.0 


NRAPR 


-297.7 


143.5 


57.0 


-28.0 


912.2 


680.1 


215.7 


232.2 


Ska25s20 


-240.9 


131.2 


2.6 


0.5 


42.3 


16.9 


-33.8 


25.4 


Ska35s20 


-240.5 


135.4 


0.0 


0.0 


-0.12 


-0.04 


-0.1 


-0.07 


SKRA 


-283.8 


152.5 


43.7 


-33.1 


699.7 


615.0 


445.5 


84.7 


SkTl 


-237.1 


122.3 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


SkT2 


-236.9 


122.1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


SkT3 


-236.9 


119.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Skxs20 


-240.6 


148.9 


4.8 


-4.9 


76.9 


77.7 


79.4 


-0.8 


SQMC650 


-263.8 


168.0 


34.3 


-36.8 


548.6 


568.3 


607.8 


-19.7 


SQMC700 


-269.6 


150.6 


39.3 


-29.2 


629.5 


548.6 


386.8 


80.9 


SV-sym32 


-253.7 


181.6 


14.5 


-37.4 


232.0 


415.3 


781.8 


-183.3 
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TABLE X: Brief compilation of methods used in fitting of the Skyrme interactions consistent with 



the macroscopic 


constraints and main 


data used in the fit. For full explanation and details, see 


original papers. 






Force 


Method 


Data used for fit 


GSkI,GSkII 


simulated annealing 
method 


ground state properties of normal and exotic nuclei: 
binding energies, charge radii, 

radii for ld^y2^^d \ij ^2 neutron orbits in "'""O and ^"'"Ca, 
breathing mode energies for ^^Zr and ^'^'^Pb, 
single-particle energies in ^"^^Pb, 
EoS of pure neutron matter [209] 


KDEOvl 


simulated annealing 
method 


ground state properties of normal and exotic nuclei: 
binding energy, charge radii and spin-orbit splitting, 
radii for ldgy2 and If 7/ 2 neutron orbits in "'"'''O and '^"'"Ca, 
breathing mode energies, critical density Pcr , 
positive slope of the symmetry energy up to 3po; 
enhancement factor associated with GDR, Landau parameter 


LNS 


Brueckner-Hartree-Fock 
with 2- and 3-body forces 
(homogeneous matter) 
Hartree-Fock 
(finite nuclei) 


nucleon effective mass in SNM and ANM, 

energy per particle in SNM and ANM as function of density and 

proton neutron asymmetry; 

constraint on the Landau parameter Got 

surface properties of selected magic and semi-magic nuclei, 
spin-orbit splitting i/2"P3/2 "'"'^O 


MSLO 


MSL [1081 model 


nuclear matter properties: 
Po , Eo , Ko , m* y , m* ^ , 

J, L, Gq, Gv [101]: surface symmetry energy. 
Landau parameters Go(Po)i Gi(po), G^{po), G'j^(po) 


NRAPR 


APR and Skyrme EoS 
of NM (a comparison) 


density dependence of effective masses as predicted by APR, 

spin-orbit splitting from charge radii and binding energies of ■^'^^Pb, ^'^Zr, ^"^Ca 


Ska25s20 
Ska35s20 


Hartrce-Fock 


binding energy, charge radii and single-particle energies 
Friedman- Pandharipande EoS 


SKRA 


Brueckner-Hartrce-Fock 
with relativistic corrections 
and 3-body forces 


nuclear matter properties: 

Po, -Eo, Xo, 

ground state properties of finite nuclei 


SkTl 

SkT2 
SkT3 


Hartrce-Fock 

Extended Thomas-Fermi 


nuclear masses and radii 

liquid drop mass formula constraints 

The same as SkTl + increase spin-orbit strength 

The same as SkTl -|- change in gradient symmetry term 


Skxs20 


Hartrce-Fock 


binding energy of doubly magic and semi-magic nuclei, rms charge radii, 
single™ particle energies, binding energy difference Ni- Ca, 
effective mass constrained to be close to unity up to 10 po 


SQMC650,700 


Quark-Meson 
Coupling model 
Hartrce-Fock 


comparison for the Skyrme (SkM*) and QMC Hamiltonians; 
free parameter: mass of the a meson 650 and 700 MeV 


SV-sym32 


Hartree-Fock 


properties of finite nuclei (energies, radii and surface thickness), 
energies of giant resonances, 

systematic variations of selected nuclear matter properties 
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FIG. 1: (color online) Constraint SMS: Pressure in SNM as a function of density up to 3^ as 
predicted by consistent Skyrme parameterizations. The shaded area in the region 2 < ^ < 4.6 is 



taken from Ref. 



78| 
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FIG. 2: (color online) Constraint SM4: Pressure in symmetric nuclear matter as a function 



of density in t^ 
cussed in Ref. 



le region 1.2 < ^ < 2.2. The shaded area represents an educated guess dis - 



80( 1 around the pressure-density relation available from kaon production data 



8ii. 



Dashed line extrapolates the pressures consistent with GMR data to higher densities in the region 



1.2 < ^ < 1.7. 

Po 
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FIG. 3: (color on line) Constraint PNMl: Energy per particle in PNM as a function of density. 
The grey bands were based on the Ref. [40] for two different parameterizations (dashed red and 
dashed-dot green lines) . The region inside the solid blue line illustrates the universal limit constraint 
given by Eq. (130]) [39 ]. For more explanation see text. 
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FIG. 4: (color online) EoS of low density pure neutron matter: (a) band defined by results of 



theoretical calculations summarized in Ref. 



42l |. See the reference for explanation of the legend. 



(b) the same as in Fig. [3l but with the additional band (a) included. 
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FIG. 5: (color online). Constraint PNM2: Pressure in the PNM as a function of density as calcu- 



lated by consistent Skyrme param eteri zations up to 3^. The bands are in the region 2 < ^ < 4.6. 



For detailed explanation see Ref. 



11C|. 
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FIG. 6: (color online) Constraints on symmetry energy So its first derivative L and Po=Ppnm(Po)) 



all at saturation density, as derived from HIC 



90|, PDR [9l|,l9a], IAS [951] and FRDM [89|. Pre- 



dictions of the consistent Skyrme parameterizations lie all within the blue dashed rectangle. For 
full explanation see text and jo?]. 
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FIG. 7: (color online) Density dependence of neutron (left panel) and proton (right panel) effective 
mass in /3-equilibrium matter as calculated by Skyrme interactions consistent with the macroscopic 
constraints. For more detail see text. 
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FIG. 8: (color online) Landau parameters calculated by CSkP sets, which passed microscopic 
constraints derived from the effective mass considerations, in SNM. See text for more explanation. 
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FIG. 9: (color online) The same as Fig. [8l but for PNM. 
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FIG. 10: (color online) Density dependence of the symmetry energy 5 as a function of j- as 
calculated by Skyrme interactions consistent with macroscopic constraints. 
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FIG. 11: (color online) Energy per particle in PNM and SNM a function of particle number density 
p for three selected Skyrme parameterizations Skxs20, SQMC700 and GSkll. 
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FIG. 12: (color online) Gravitational mass vs radius for cold noi i-rot ational neutron stars as 



calculated using a Skyrme EoS aug men ted by BPS EoS at low density 



14g] (left panel) and matched 



by a FQMC EoS at high densities 149[] and by BPS EoS at low densities (right panel). The dashed 
lines indicate the limits on the maximum mass of the most massive neutron star observed up- 
to-date [147l |. Only Skyrme parameterizations which are consistent with both macroscopic and 



microscopic constraints are used. For more explanation see text. 
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FIG. 13: (color online) Relation between the gravitational mass, Mg, for the selected Skyrme 
models, and the corr espo nding baryonic mass, Mb- The boxes represent constraints derived by 



Podsiadlowski et al. 



150l | (full line box) and more recently by Kitaura et al. [ 15l[ | (dashed line 



box) based on the proposed properties of system J0737-3039, as discussed in the text. Results are 
shown for Skyrme interaction consistent with both, macroscopic and microscopic constraints. 
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FIG. 14: (color online) Density dependence of the energy per particle, resulting only from the 
contributions of the non-standard terms in Eq. ([T]) for all non-standard Skyrme parameterizations 
used in this work. 



86 



